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EXECUTIVE  SUMMARY 


I.  Summary  of  the  Program 

This  report  summarizes  the  research  that  was  carried  out  under  the  Joint  Services  Electronics 
Program  at  the  University  of  Illinois  at  L'rbana-Champaign  during  the  period  April  1.  1986. 
through  March  31.  1987  Since  the  new  3-vear  JSEP  contract  began  on  October  1.  1986.  this 
reporting  period  covers  the  last  six  months  of  the  previous  JSEP  contract,  as  well  as  the  first  six 
months  ol  the  present  one.  The  JSEP  contract  that  began  on  October  1.  1986.  contains  22  work 
units  (Unit  3  ot  the  new  contract  was  not  funded),  which  is  a  reduction  I  rom  the  26  work  units  in 
t he  previous  contract.  Two  previously  supported  units  in  the  surlace  physics  area  were 
dropped:  (1)  the  unit  on  acoustic  charge  transport  (Professor  Hunsinger)  is  now  funded  by 
DARPA.  and  (2)  Professor  Davidson  has  become  Associate  Director  of  the  Center  for  Supercom¬ 
puter  Research  and  Development.  His  work  on  computer  architecture  is  now  funded  by  DOE.  \SF. 
and  private  industry.  The  old  unit  (Unit  3)  on  vapor  phase  growth  and  characterization  was 
replaced  by  the  new  Unit  5  on  MOCVD-grown  heterostructure  electronic  devices  under  the  direc¬ 
tion  of  Professor  J.  J.  Coleman.  In  the  new  contract.  Professor  Thornton  has  been  added  to  inves¬ 
tigate  chemistry  processes  in  plasma  discharges  (Unit  9).  and  Professor  I. v ding,  a  new  recruit,  was 
added  to  Unit  8  (formerly  Unit  10)  to  perform  some  key  experiments  supporting  the  charge  den¬ 
sity  wave  theory  work  under  the  direction  of  Professor  Tucker.  The  prev  ious  unit  on  hierarchical 
design  aids  was  replaced  by  the  new  Unit  7  on  high-performance  integrated  circuits.  Two  new 
faculty.  Professors  Kang  and  Rao.  have  replaced  Professor  Trick,  who  has  become  the  new  Head  of 
Electrical  and  Computer  Engineering  and  another  faculty  member  who  is  no  longer  with  the 
university  (Professor  Mayeda). 

In  the  electromagnetics  area.  Professors  Chuang  and  Lo  are  new-  to  the  unit  on  monolithic  mil¬ 
limeter  wave  ICs  (Unit  11).  Professor  Chuang  is  a  new  faculty  member  from  MIT  with  a  unique 
background  in  both  semiconductors  and  electromagnetics.  In  the  new  contract.  Unit  12  on  inverse 
scattering  under  the  direction  of  Professor  Mittra  was  expanded  due  to  increased  interest  in  this 
subject. 

In  the  information  systems  area,  three  new  faculty  were  recruited  in  fault-tolerant  comput¬ 
ing  (Professors  Fuchs.  Banerjee.  and  Iyer),  and  a  new  Unit  17  on  parallel  VLSI  structures  for  sen¬ 
sor  array  processing  has  been  added  with  two  new  faculty.  Professors  Arun  and  Wah.  Finally,  a 
total  of  six  new  faculty  in  the  areas  of  adaptive  and  nonlinear  systems,  communication  systems, 
and  statistical  signal  processing  were  added  at  the  beginning  of  the  new  contract  in  October  1986 
(Professors  Kumar.  Medanic.  Grizzle.  Poolla.  Arikan.  and  Barron). 

Professor  Jenkins  has  served  as  Acting  Director  of  the  Coordinated  Science  Laboratory  during 
the  last  academic  year  to  replace  Professor  Trick,  and  both  Jenkins  and  Trick  have  been  sharing  the 
responsibilities  of  JSEP  principal  investigator  since  August  1,  1986.  Recently.  Professor  Jenkins  has 
been  named  as  the  new  Director  of  CSL  and  will  assume  the  position  on  a  permanent  basis  on  May 
21.  1987. 

II.  Most  Outstanding  Accomplishments 

There  are  three  particular  items  that  we  would  like  to  cite  as  the  most  significant  accomplish¬ 
ments  on  the  JSEP  Program  during  this  reporting  period,  each  of  which  is  described  briefly  below. 
These  are  the  development  of  a  new  heterostructure  hot  electron  diode,  new  error-control  coding 
and  spread-spectrum  communication  results  that  have  been  incorporated  into  the  Army's 
SINCGARS  radio,  and  the  successful  design  and  fabrication  of  a  VLSI  integrated  circuit  that  con¬ 
tains  a  programmable  digital  filter  module  whose  design  is  based  on  modern  residue  number  arith¬ 
metic  techniques. 


The  Heterostructure  Hot  Electron  Diode 

Professors  J.  J.  Coleman  and  K..  Hess,  together  with  graduate  students  T.  K.  Higman  and 
\I  A  Kmanuel  have  discovered  a  new  switching  mechanism  in  semiconductor  heterostructures 
which  has  resulted  in  the  development  of  a  new  electronic  device,  the  heterostructure  hot  electron 
diode  (  H2ED).  The  operation  of  the  H2ED  is  based  on  a  transition  between  two  current  conduction 
modes  in  a  two-terminal  structure,  a  low  current  tunneling  mode,  and  a  high  current  thermionic 
emission  mode.  Sw  itching  between  these  modes  gives  rise  to  an  S-shaped  negative  differential  resis¬ 
tance  in  the  current  versus  voltage  characteristic.  These  devices  have  application  as  a  microwave 
oscillator.  Microwave  characterization  is  presently  in  progress,  and  preliminary  results  have 
demonstrated  test-fixture-limited  gain  to  IS  GHz.  If  no  parasitic  effects  are  encountered,  and  if  the 
oscillation  frequency  is.  in  fact,  electron  transit  time  limited,  oscillation  in  the  HXX)  GHz  range 
may  be  possible. 

This  work  represents  a  synergistic  relationship  between  Professor  Hess'  theoretical  work  sup¬ 
ported  under  Unit  2  and  Professor  Coleman's  experimental  work  supported  under  Unit  5.  Profes¬ 
sor  Hess  developed  the  theoretical  model  of  the  new  device  and  predicted  its  characteristics.  Pro¬ 
fessor  Coleman  built  the  new  device  in  his  laboratory,  measured  the  physical  characteristics  of  his 
experimental  sample,  and  is  presently  continuing  its  further  characterization  and  development.  It 
was  announced  in  a  press  release  by  the  University  of  Illinois  at  Urbana-Champaign  to  the  UP1 
w  ire  service  that  a  patent  disclosure  has  been  filed  on  this  new  device,  and  a  number  of  publications 
have  either  appeared  or  will  be  forthcoming  soon  describing  the  technical  details  of  the  H2ED  (see 
Unit  5.  Research  Summary). 

Error-Control  Coding  and  Spread-Spectrum  Communication  Techniques  for  the  Army’s 
SINCGARS  Radio 

Several  results  of  CSL  research  in  error-control  coding  and  spread-spectrum  communications 
have  recently  been  incorporated  into  the  SINCGARS  radio  (Single  Channel  Ground/Airborne  Radio 
Subsystem).  This  radio,  the  Army's  newest  tactical  radio,  was  developed  by  ITT 
Aerospace  Optical  Division  for  voice  and  data  communications.  It  operates  in  both  a  single-channel 
mode  and  a  frequency-hop  spread-spectrum  mode.  In  addition,  an  overlay  is  being  developed  by 
ITT  and  SRI  International  to  enable  SINCGARS  radios  to  operate  as  packet  radios  in  a  mobile  com¬ 
munication  network. 

Research  results  obtained  by  CSL  Professors  M.  B.  Pursley  and  D.  V.  Sarwate  and  their  gra¬ 
duate  research  assistants  were  employed  by  Professor  Pursley  to  design  the  synchronization 
sequence  that  is  utilized  for  timing  acquisition  in  the  SINCGARS  receiver.  In  addition,  for  the 
packet  radio  version  of  SINCGARS  (known  as  the  SINCGARS  Packet  Overlay).  Professor  Pursley 
developed  the  error-control  method  which  combines  Reed-Solomon  codes  with  side  information 
derived  from  test  patterns.  The  side  information  is  used  to  erase  unreliable  data  at  the  output  of 
the  demodulator  so  that  the  input  to  the  decoder  is  a  combination  of  data  symbols  and  special  eras¬ 
ure  symbols.  The  code  enables  the  receiver  to  correct  errors  and  erasures  caused  by  noise  and 
interference  in  the  communication  channel.  The  combination  of  coding  and  side  information  pro¬ 
vides  the  capability  for  a  large  number  of  SINCGARS  radios  to  operate  simultaneously  in  the  same 
frequency  band,  and  it  also  improves  the  radio’s  ability  to  communicate  in  the  presence  of  various 
types  of  electromagnetic  interference. 

The  decoding  algorithm,  which  was  developed  and  implemented  in  software  by  Professor 
Sarwate.  is  basically  a  mathematical  description  of  the  computations  that  are  required  to  correct 
both  errors  and  erasures.  The  actual  computations  prescribed  by  the  algorithm  are  carried  out  in  a 
microprocessor  contained  in  the  radio.  The  resulting  decoder  is  capable  of  handling  data  rates  in 
excess  of  16.000  data  symbols  per  second,  the  maximum  signaling  rate  of  the  SINCGARS  radio. 


Professor  Pursley  also  provided  the  analytical  tools  to  evaluate  the  performance  of  the  error- 
control  scheme  and  assess  its  impact  on  the  overall  communication  network  performance.  Analyti¬ 
cal  results  from  this  effort,  and  other  investigations  conducted  by  Professor  Pursley  and 
Mr  J  R  McChesney  ol  ITT.  were  used  subsequently  in  the  SINCGARS  packet  radio  network 
simulation  developed  for  the  Army  by  SRI. 

The  research  leading  to  the  development  of  the  approach  used  for  error  control  in  the 
SINCGARS  radio  was  funded  bv  JSFP  during  the  early  stages,  funded  by  ARO  during  later  stages, 
and  developed  into  the  actual  SINCGARS  application  under  LI  contracts  with  ITT  and  individual 
consulting  arrangements  with  SRI  and  ITT. 

Design  and  Fabrication  of  a  VLSI  Digital  Filter  Module 

For  a  number  of  years  both  \SF  and  JSFP  have  supported  theoretic  studies  at  the  Coordinated 
Science  Laboratory  on  the  use  of  residue  number  arithmetic  (R.\S)  for  the  design  of  high-speed 
VI  SI  circuits  for  the  digital  signal  processing  functions  that  are  required  in  modem  communication 
and  radar  systems.  During  the  last  year.  Professors  Jenkins.  Arun.  and  Wah  have  carried  out  the 
successful  design  and  fabrication  of  a  semicustom  VLSI  digital  filter  module  which  demonstrates 
both  the  feasibility  and  desirability  of  using  residue  number  theory  design  techniques  in  practice 
(Unit  17).  A  semicustom  integrated  circuit  module  was  designed  for  an  RMS  digital  filter  using  the 
IBM  MVISA  C.AD  system  that  is  located  on  the  campus  of  the  University  of  Illinois  at  Urbana- 
Champaign.  A  programmable  4-bit  module  was  designed  as  a  basic  building  block  for  a  finite 
impulse  response  digital  filter  that  has  order  8  and  14  bits  of  arithmetic  precision.  The  module  was 
designed  using  finite  field  logarithm  additions  to  completely  eliminate  conventional  multipliers  on 
the  chip  and  thereby  improve  speed  and  circuit  density.  The  chip  design  was  done  with  MVISA 
using  a  standard  cell  approach,  and  MVISA  was  then  used  to  simulate  the  hardware  operation  and 
to  estimate  physical  parameters.  The  chip  was  subsequently  fabricated  at  the  IBM  Manassas  facil¬ 
ity.  and  10  samples  were  recently  returned  to  the  University  of  Illinois  at  Urbana-Champaign  for 
testing  and  characterization.  MVISA  estimates  that  the  4.2  mm  chip  consumes  89  mW  and 
operates  at  a  system  cycle  frequency  of  10  mHZ,  which  corresponds  to  a  data-cycle  frequency  of 
1.2  mHz  The  design  used  529  out  of  a  possible  560  available  standard  cells.  Correct  operation  of 
the  circuits  has  not  yet  been  verified  since  they  were  received  only  two  weeks  ago.  However,  since 
the  circuits  were  designed  with  testability  in  mind,  and  due  to  the  modular  nature  of  the  RNS 
architecture,  it  is  believed  that  testing  will  be  a  straightforward  matter. 

The  development  of  this  VLSI  digital  filter  module  would  not  have  been  possible  without  the 
theoretical  support  that  was  provided  by  NSF,  VHSIC.  and  JSEP  for  many  years,  the  synergistic 
research  interaction  between  the  Computer  Systems  Group  (Units  13  and  14)  and  the  Digital  Signal 
and  Image  Processing  Group  (Units  16  and  17).  and  the  support  of  private  industry  (IBM)  who 
served  as  a  silicon  foundry  for  chip  fabrication.  This  work  was  recently  presented  at  the  1987 
International  Symposium  on  Circuits  and  Systems  and  will  be  described  in  further  detail  in  several 
forthcoming  publications  (see  Unit  17).  It  is  expected  that  the  interdisciplinary  team  effort  that 
was  put  together  for  this  project  will  serve  as  a  model  for  future  research  projects  of  this  type. 
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WORK  UNIT  NUMBER  1 


TITLE:  Crystal  Growth  from  the  Vapor  Phase  and  Controlled  Doping  of  Equilibrium  and  Meta¬ 
stable  Semiconductor  Alloys:  Ion-Surface  Interactions 


SENIOR  PRINCIPAL  INVESTIGATORS: 

J  I  Greene.  Research  Professor 
S.l  Shah.  Research  Associate 

SCIENTIFIC  PERSONNEL  AND  TITLES: 

P.  Fons.  Research  Assistant 
S.  Gorbatkin.  Research  Assistant 
B.  Kramer.  Research  Assistant 
D.  Lubben.  Research  Assistant 
D  McIntyre.  Research  Assistant 
D.  Mei.  Research  Assistant 
M.A.  Ray.  Research  Assistant 
A.  Rockett.  Research  Assistant 
L.  Romano.  Research  Assistant 


SCIENTIFIC  OBJECTIVE: 

The  primary  objective  of  this  research  program  is  to  investigate  energetic  particle-surface 
interactions  which  control  the  nucleation  and  growth  kinetics,  chemistry,  and  physical  properties 
of  alloy  semiconductors  during  vapor  phase  crystal  growth  by  UHV  ion  beam  sputtering  and 
accelerated-beam  molecular  beam  epitaxy.  In  both  of  these  growth  techniques,  low  energy  ion- 
surface  interactions  allow  an  efficient  coupling  of  kinetic  energy  to  the  growth  surface  upon  con¬ 
densation.  thereby  altering  the  surface  reactivity  as  well  as  adsorption  and  adatom  diffusion  kinet¬ 
ics  allowing  single  crystal  film  growth  at  lower  temperatures,  much  more  precise  control  over 
dopant  incorporation  probabilities  and  depth  distributions,  and  the  growth  of  unique  metastable 
alloys.  This  work  is  being  pursued  from  both  an  analytical  and  an  experimental  point  of  view  in 
order  to  establish  a  detailed  understanding  of  the  fundamental  film  growth  mechanisms. 

SUMMARY  OF  RESEARCH: 

Surface  Reaction  Mechanisms,  Elemental  Incorporation  Probabilities,  and  Depth  Distribu¬ 
tions  of  Accelerated  and  Thermal  Dopants  in  Semiconductor  Films  Grown  from  the  Vapor 
Phase 

As  reported  last  year,  we  have  developed  a  general  time-dependent  model,  which  combines 
thermodynamic  and  kinetic  elements,  for  describing  the  incorporation  of  dopants  into  films  during 
deposition.  The  model  accounts  for  dopant  surface  segregation  and  allows  elemental  incorporation 
probabilities  <J  and  depth-dependent  concentration  profiles  C(x.t)  to  be  calculated  as  a  function  of 
experimental  parameters  such  as  film  and  dopant  material,  deposition  rate,  incident  dopant  flux, 
growth  temperature,  etc. 

Calculated  profiles  from  arbitrarily  complex  doping  schedules  for  both  thermal  and 
accelerated  dopants  are  in  excellent  agreement  with  experimental  results  for  a  variety  of  dopants  in 
MBE  Si.  MBE  and  sputter-deposited  GaAs,  and  MBE  Ga|_xAlxAs.  Moreover,  we  have  used  the 
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model  to  predict  critical  temperatures  for  transitions  in  dopant  segregation  regimes,  dopant-induced 
surface  structural  phase  transitions,  and  changes  in  dopant/surface  reaction  paths  leading  to  large 
modifications  in  surface  binding  energies  and  hence  incorporation  probabilities.  All  of  these  effects 
have  now  been  observed  experimentally  Additional  terms  have  been  added  to  the  model  to 
account  for  trapping  (i.e.  low-energy  implantation)  and  enhanced  segregation  associated  with  the 
use  of  accelerated  dopants.  The  extended  model  has  been  shown  to  work  very  well  for  acceleration 
energies  of  a  few  hundred  eV  or  more  where  low-energy  implantation  is  the  dominate  effect. 

Selected  results  from  our  experimental  investigations  are  briefly  summarized  below. 

•  I  he  Si  incorporation  probability  <xs.  and  segregation-induced  profile  broadening  SSl  of  nomi¬ 
nally.  abrupt  Si  doping  modulations  during  the  growth  of  MBE  Ga|_xAlyAs  have  been  determined 
as  a  1  unction  ol  IL.  x.  and  As2  overpressure  PAs  during  growth.  Si  segregation  in  Ga|_xAlxAs.  while 
quite  small  compared  to  other  dopants  such  as  Sn.  was  (contrary  to  popular  belief)  measurable.  8Sl 
ranged  form  —  7  nm  at  600°C  to  25  nm  at  72 5°C  and  segregation  was  kinetically  limited  over  this 
Ts  range  Ss,  increased  slightly  with  x,  but  varying  PAs  by  a  factor  of  three  had  no  measurable 
effect.  A  lower  limit  for  the  Gibbs  free  energy  of  segregation  was  found  to  be  0.4  eV.  cr^  varied 
from  1  at  Ts  ^  600 °C  to  ~~  0.5  at  725°C  and  was  not  a  strong  function  of  either  x  or  PAs.  The  sur¬ 
face  binding  energy  of  Si  was  estimated  to  be  3  ±  0.1  eV. 

•  The  incorporation  probability  of  thermal  In.  a  deep  acceptor,  in  \1BE  Si(  100)  was  undetect¬ 
able  at  normal  growth  temperatures  and  was  of  the  order  of  10~s  at  Ts  =  475°C.  This  low  incor¬ 
poration  probability  results  from  the  combined  effects  of  a  strong  tendency  for  In  atoms  to  follow 
the  advancing  Si  film  surface  and  the  relatively  high  desorption  rate  of  In  at  Si  growth  tempera¬ 
tures.  The  steady-state  In  surface  coverage  0|„  on  the  growing  Si  surface  accumulates  to  cause  the 
surface  reconstruction  to  change  from  the  initial  Si(2xl )  to  In-stabilized  (3x4),  corresponding  to  0,n 
>  0.1  monolayer  (ML),  and  then  to  In(3x4)  with  SiC 310)  facets  at  0In  >  0.7  ML  without  any 
detectable  In  incorporation  in  the  Si  film. 

Lsing  a  single-grid.  UHV  compatible  ion  source  of  our  own  design  to  provide  accelerated  ln  + 
doping  beams,  we  have  obtained  crIn  values  of  0.02-0.7  in  Si(100)  films  grown  at  800°C  with  In+ 
acceleration  energies  E|n  between  50  and  400  eV.  Combined  SIMS  and  T-dependent  Hall  and  resis¬ 
tivity  measurements  carried  out  on  films  grown  with  E!n  =  200  eV  showed  that  In  was  incorporated 
substitutionallv  into  electricallv-active  sites  over  the  entire  concentration  range  examined. 
10lh— 101'’  cm-3  (well  above  equilibrium  solid-solubililv  limits!).  The  acceptor  level  ionization 
energy  was  165  meV.  in  agreement  with  bulk  Si  results,  and  there  was  no  evidence  for  the  111 
meV  level  associated  with  In-C  complexes  that  is  reported  for  annealed  In-implanted  Si.  Carrier 
mobilities  were  higher  than  the  best  results  reported  for  In-implanted  Si  and  above  Irvin-curve 
values  for  bulk  Si.  There  was  no  evidence  for  residual  lattice  damage  in  the  films  until  E|n  was 
increased  to  500  eV  while  Ts  was  decreased  to  600 “C  where  incomplete  In  substitutionality  and  a 
loss  of  ~  30%  in  mobility  were  observed. 

•  We  previously  determined  and  modeled  the  detailed  mechanism  and  the  kinetics  of  the  dis¬ 
sociative  chemisorption  of  Sb4  on  SiC  100)2x1  surfaces  (see  earlier  reports  and  Barnett.  Winters,  and 
Greene.  Surf.  Sci  165.  303  (1986))  and  calculated  thermal  Sb  incorporation  behavior  in  MBE  Si 
(see  Barnett  and  Greene.  Surf.  Sci.  151.  67  (1985)).  We  have  recently  begun  experiments  with  a 
new  ion  gun  to  investigate  the  incorporation  of  low-energy  accelerated  Sb+  ions  during  MBE  Si. 

Low-Energy  Ion/Surface  Interactions  During  Crystal  Growth 

Low-energy  ion-irradiation  of  the  film  during  deposition  can  be  used  to  provide  an  efficient 
coupling  of  kinetic  energy  to  the  growth  surface  in  order  to  alter  surface  reactivity  as  well  as 
adsorption,  adatom  diffusion,  nucleation.  and  growth  kinetics.  Such  effects  have  been  reported  by 
laboratories  across  the  world  using  a  variety  of  film  growth  techniques  including  accelerated-beam 
MBE.  Pa-CVD,  and  sputter  deposition  (see  Greene.  Solid  State  Technology  April.  1987  for  a  recent 
review).  This  in  turn  can  lead  to  lower  epitaxial  temperatures,  more  control  over  elemental  incor¬ 
poration  probabilities,  and  the  growth  of  unique  metastable  alloys.  Our  group  has  been  modeling 


ion/surface  interaction  effects  which  are  common  to  all  of  these  techniques  and  have  carried  out  the 
first  definitive  experiments  under  well-controlled  conditions  to  probe  fundamental  mechanisms. 
Some  highlights  are  described  briefly  below. 

•  We.  and  others,  have  predicted  that  one  of  the  important  roles  of  low-energy  ion  irradiation 
during  the  initial  stages  of  film  nucleation  is  to  provide  preferred  adsorption  sites  (i.e.  non-random 
nucleation  sites).  We  have  recently  developed  experimental  techniques  using  thermally  stimulated 
desorption  (TSD)  and  modulated-beam  mass  spectrometry  (MBMS)  to  actually  measure  such 
effects  for  the  first  time.  Si(  100)2x1  surfaces  at  Ts  =  425°C  were  irradiated  in  separate  experi¬ 
ments.  during  and  after  deposition  of  1  Ml.  of  Sb.  In  both  cases,  in  addition  to  the  normal  thermal 
Sb  binding  energy  site  at  2.33  eY  (which  we  had  reported  previously),  we  observed  an  additional 
higher  binding  energy  site  at  2.6  eV  The  Sb  coverage  in  the  E2  site  increased  w  ith  ion  dose  and  a 
model  was  developed  that  provided  a  good  fit  to  the  data  and  allowed  us  to  calculate  Sb  trapping 
probabilities. 

•  We  have,  for  the  first  time,  investigated  the  nucleation  and  growth  of  films  deposited  from 
accelerated  beams  under  well  controlled  conditions  in  L'HV.  Transmission  electron  microscopy  was 
used  to  directly  observe  the  kinetics  of  film  growth  from  partially  ionized  In  beams  on  electron- 
transparent  Si3\4  windows.  The  use  of  acceleration  energies  ranging  from  0  to  300  eV  were  found 
to  dramatically  alter  nucleation  kinetics  leading  to  a  progressive  shift  toward  larger  average  island 
diameters  with  a  more  uniform  size  distribution.  Secondary  island  formation  was  almost  entirely 
suppressed  due  primarily  to  ion-irradiation-induced  dissociation  of  small  clusters  and  enhanced 
adatom  diffusion.  However,  we  believe  that  the  mechanism  for  enhanced  adatom  diffusion  is  quite 
different  than  the  one  that  is  popularly  held. 

Crystal  Growth  and  Physical  Properties  of  Single-Crystal  Metastable  Semiconductors 

We  have  carried  out  the  first  detailed  studies  of  the  growth  and  physical  properties  of  new 
single-crystal  melastable  semiconductors.  The  key  feature  stabilizing  the  growth  of  these  material 
is  the  controlled  use  of  low-energy  ion  bombardment  to  modify  the  incorporation  probabilities  of 
the  matrix  elements  and  to  promote  dynamic  collisional  mixing  during  deposition. 

•  Single-phase  metastable  Gei_xSnx  alloys  with  a  diamond-cubic  structure  have  been  grown 
with  Sn  concentrations  up  to  15  at  %  (the  maximum  equilibrium  solid  solubility  is  —  1  at  %). 
Films  with  x  ^  0.08  grown  on  Ge(100)  and  GaAs(lOO)  substrates  were  single  crystals  while  films 
grown  on  amorphous  glass  were  polycrystalline  with  a  (220)  preferred  orientation.  The  allowable 
range  in  obtaining  single-phase  crystalline  alloys  was  a  function  of  x  (the  thermodynamic  driving 
force  for  phase  separation  increases  with  increasing  Sn  concentration)  and  the  ion  flux  and  accelera¬ 
tion  energy.  A  growth  phase  map  plotted  as  a  function  of  x.  Ts,  and  E  was  determined.  Raman 
spectra  of  Gei_xSnx  alloys  exhibit  a  two  peak  density  of  states  which  can  be  explained  by  assuming 
the  Sn  substitution  introduces  disorder  in  the  bonding  direction  and  removes  the  degeneracy  of  the 
opticr  '  branches. 

•  GaAsj_xSbx  alloys,  of  interest  for  fabrication  of  sources  and  detectors  in  the  ^  1  fim  range 
for  fiber-optics  communications,  have  been  grown  on  GaAs(  100)  substrates  at  compositions  across 
the  entire  pseudobinary  range.  Under  equilibrium  conditions,  the  alloy  system  exhibits  a  miscibil¬ 
ity  gap  from  x  —  0.2  to  0.9  at  the  growth  temperature.  X-ray  diffraction,  optical  absorption,  and 
Raman  spectra  showed  that  our  samples  were  single-crystal,  single-phase,  alloys  for  all  x.  Alloy 
lattice  constants  varied  linearly  with  concentration  while  the  direct  T-point  optical  bandgap 
showed  significant  bowing  with  a  minimum  near  x  =  0.8.  Raman  spectra  exhibited  a  two-mode 
behavior  throughout  the  alloy  system  including  a  local  mode  of  As  in  GaSb  and  a  resonant  mode  of 
Sb  in  GaAs.  This  is  contrary  to  previous  reports  of  a  “one-two"  mode  behavior  in  OMCVD 
GaAs]_xSbx  for  which  spinoidal  decomposition  was  suspected.  Moreover,  the  Raman  peaks  from 
our  films  had  peak  broadenings  and  assvmmetries  much  less  than  these  observed  in  the  OMCVD 
samples  and.  in  fact,  less  than  those  of  VIBE  Ga1_xAlxAs. 


7 


y 

[9] 

v 

H 

[10] 

y 

y 

[11] 

V, 

[12] 

‘  s- 

V 

i 

[13] 

r. 

r 

JSEP-SPONSORED  PUBLICATIONS 

S  A.  Barnett,  H  F  Winters,  and  J  F  Greene  "The  interaction  of  Sb4  molecular  beams  with 
Sit  100)  surfaces:  Modulated-beam  mass  spectrometry  and  thermally  stimulated  desorption 
studies." Surf  Sci ..  col.  165,  p.  303,  19X6.  (JSFP  SRC) 

J  knall.  J.-F.  Sundgren.  G  V  Hansson.  and  J  F  Greene,  "in  overlayers  on  clean  Sit  100)2x1 
Surface  structure,  nucleation.  and  growth."  Surf  Sci  .  col.  166.  p.  512.  19X6  (JSFP  SRC) 

A  Rockett.  J  knall.  \1  A.  Hassan.  J.-h.  Sundgren.  S  A  Barnett  and  J  F  Greeene.  "Ther¬ 
mal  anil  accelerated  (  $  200  e\  )  In  iloping  ot  Sit  100)  lavers  during  molecular  beam  epitaxy." 
J  Vac  Sci  Technol  A.  col.  4.  P.  900.  19X6.  (JSFP  SRC) 

A.  Rockett.  J  klem.  S.  A.  Barnett.  J.  E  Greene,  and  H.  Vlorkoc.  "Si  incorporation  and  segre¬ 
gation  on  Gaj_xAlxAs  ( 100)  films  grown  bv  molecular  beam  epitaxv,"  J  Vac.  .Set  Technol  R, 
col. -4.  p.  519.  19X6.  (JSFP) 

A.  Rockett.  J.  klem,  S.  A.  Barnett.  J.  E.  Greene,  and  H  Morkoc.  “Si  incorporation  probabili¬ 
ties  and  depth  distributions  in  Ga,_xAlxAs  films  grow  in  bv  molecular  beam  epitaxy,"  J  Appl 
Phys  .  col  59.  P  2777.  19X6.  (JSFP) 

T  C  McGlinn.  T  \  krabach.  \1  V  klein.  G  Bajor.  J  F  Greene.  B.  kramer.  S.  A  Barnett. 
A  Lastras.  and  S  Gorbatkin.  "Raman  scattering  and  optical  absorption  studies  of  the  meta- 
stable  alloc  svstem  Ga.As,_xSbx." /Vnr  Rev  R.vol  33.  p.  X396.  19X6.  (JSF.P) 

f 

J  F  Greene.  T  Motooka.  J.-E.  Sundgren.  A.  Rockett.  S.  Gorbatkin.  D.  Lubben.  and  S.  A.  Bar¬ 
nett.  "A  review  ot  the  present  understanding  of  the  role  of  ion- surface  interactions  and 
photo-induced  reactions  during  vapor  phase  crvstal  growth."  J  Crvst.  Growth,  vol  79.  p.  19. 
19X6  (JSFP  DOF  SRC) 

M  A  Hasan,  J  knall.  S  A  Barnett,  J.-E.  Sundgren.  A  Rockett  and  J.  F.  Greene. 
"Accelerated-ion  beam  doping  during  Si  growth  by  molecular  beam  epitaxy  and  Ion-enhanced 
In  film  deposition  using  low  energy  (40  to  300  eV)  In  +  Ion  source."  Vacuum,  in  press. 
(JSEP  SRC) 

S.  A.  Barnett.  H.  F  Winters,  and  J.  E.  Greene,  "Influence  of  Ion  bombardment  on  the  interac¬ 
tion  of  Sb  with  the  Si(  100)  surface,"  Surf  Sci  .  in  press.  (JSFP  SRC) 


tastable  Ge|_xSnx  allovs  on  Ge(  100)  and  GaAs(  100)  substrates."  J  Crvstal  Growth,  in  press. 
(JSEP  DOE) 


stable  semiconductors."  J  Vac  Sci  Technol..  in  press.  (JSEP  DOE) 

J  h.  Greene.  T  Motooka.  J.-E.  Sundgren,  D.  Lubben.  S.  Gorbatkin,  and  S.  A.  Barnett.  "The 
role  of  Ion/ surface  interactions  and  photo-induced  reactions  during  film  growth  from  the  va¬ 
por  phase."  /  \ucl  Instr  Methods  B.  in  press.  ( JSEP/DOE,  SRC) 


Effects  on  microchemistry  and  microstructure."  Sol  St  Technol  ,  in  press.  (JSEP  SRC-  DOE) 


[14]  D.  H  Rich.  A  Samsavar.  T.  Miller  H.  F  Lin.  T.-C  Chiang.  J.-F,  Sundgren  and  J  1  Creene 
Coordination  determination  ot  In  on  Si(  100)  from  synchrotron  photoemission  studies."  Ph.s 
Rev.  Lea.,  in  press  (JSFP  SRC) 

[15]  M.-A.  Hasan.  S  A  Barnett.  J.-F.  Sundgren.  and  )  F  Creene.  "Xucleation  and  initial  growth 
ol  In  deposited  on  Si  ,\4  using  low -energy  300  eV  )  accelerated  beams  in  ultra  high  V  acu¬ 
um.  J  \  ac  Sci  Tecimol  in  press.  (JSFP  SRC) 


9 


I 


WORK  UNIT  NUMBER  2 


TITLE:  Studies  of  Transport  Phenomena  in  Semiconductors 


SENIOR  PRINCIPAL  INVESTIGATORS: 


k  Hess.  Research  Professor 

.1  P  I  eburton  Research  Assistant  Professor 


SCIENTIFIC  PERSONNEL  AND  TITLES: 


D.  Bailey,  Research  Assistant 
J  Higman.  Research  Assistant 
S  Manion.  Research  Assistant 


SCIENTIFIC  OBJECTIVE: 

This  research  involves  the  study  of  basic  properties  of  semiconductors,  semiconductor- 
heterolayers.  new  device  concepts  and  device  simulation  Both  theoretical  and  experimental 
methods  are  employed  in  each  of  these  categories.  We  are  examining  a  variety  of  hot  electron 
phenomena  and  their  effects  on  present  and  future  device  performance,  especially  in  connection 
with  modulation  doping  The  experiment  studies  concern  mainly  electronic  transport  in  hetero- 
layers  in  high  electric  and  high  magnetic  fields. 

The  theoretical  studies  include  Monte  Carlo  simulations  of  electronic  transport  and  the 
development  of  two-dimensional  models. 


SUMMARY  OF  RESEARCH: 

Novel  Heterolayer  Devices  and  Transport  Effects 

We  have  proposed  a  new  switching  mechanism  in  a  two-terminal  semiconductor  heterolayer 
structure  which  capitalizes  on  nonlinear  electron  temperature  effects  in  adjacent  heterolayers  [l]. 
We  have  been  able  to  demonstrate  experimentally  that  this  mechanism  can  be  used  to  fabricate 
oscillators  with  theoretical  frequency  limits  in  the  Teraherz  range.  Oscillations  have  been  observed 
up  to  20  GHz  limited  by  the  experimental  set  up. 

We  have  also  investigated  dual  channel  high  electron  mobility  transistors  and  their  ultimate 
switching  speed  by  ensemble  Monte  Carlo  simulations  and  have  found  switching  speeds  superior  to 
conventional  high  electron  mobility  transistors  [3.9]. 

Furthermore  we  have  investigated  tunneling-assisted  impact  ionization  across  the 
conduction-band-edge  discontinuity  of  quantum-well  heterostructures  and  applied  the  results  to  a 
new  superlattice  structure.  We  have  considered  muhiquantum-well  structures  where  the 
quantum-well  regions  are  heavily  doped  and  the  undoped  barrier  regions  are  essentially  insulating 
Incident  hot  electrons  due  to  the  applied  electric  field  perpendicular  to  the  heterointerface  interact 
with  the  two-dimensional  electrons  confined  to  the  quantum  wells  through  the  Coulomb  force. 
The  resultant  electrons  can  either  have  enough  energy  to  get  out  of  the  wells  or  tunnel  through  the 
triangular  barriers.  A  new  analytical  approximation  for  the  impact  ionization  rate  is  given  which 


compares  favorably  with  the  numerical  results.  The  tunneling-assisted  impact  ionization  rates  and 
the  ionization  coefficients  are  also  calculated.  It  is  shown  that  the  tunneling  effect  reduces  the  ioni¬ 
zation  thresholds  and  enhances  the  ionization  rates  significantly.  This  effect  is  a  candidate  for  new 
forms  of  solid  state  devices  analog  to  photomultipliers  [f>] 

Electronic  Transport  at  High  Energies  and  Transient  Phenomena; 

We  have  continued  to  investigate  the  basic  mechanisms  responsible  for  the  conductance  oscil¬ 
lations  in  the  "Hickmolt"  experiment.  We  have  been  able  to  complete  a  model  based  on  phononiza- 
tion  ol  neutral  donors  in  the  active  GaAs  layer  ol  the  semiconductor  structure.  Recent  experiments 
using  far  mtrared  spectroscopy  have  demonstrated  the  occurrence  of  ionization  in  the  oscillation 
cycles,  which  definitively  excludes  interpretations  based  on  “exotic"  effects  such  as  magnetopolaron 
or  Ihm  s  model.  Impact  ionization,  as  proposed  bv  haves  and  coworkers,  has  been  claimed  to  occur 
via  a  magnetic-impurity  resonance  (MIR)  effect.  We  have  demonstrated  that  the  interpretation  is 
questionable.  Both  phonon  and  impact-ionization  mechanisms  seem  so  far  undistinguishable. 
Further  experiments  are  necessary  to  identify  the  fundamental  process  causing  the  current  oscilla¬ 
tions.  [4.5]. 

A  second  investigation  (experimental)  of  low-temperature,  high-field  transport  has  been  deal¬ 
ing  with  the  energy  loss  rates  of  electrons  in  modulation-doped  AlGaAs  GaAs  heterojunctions  at 
low  temperatures.  At  the  temperatures  of  this  experiment,  the  energy  loss  is  due  to  acoustic  pho¬ 
non  scattering,  mostly  through  the  deformation  potential  interaction.  We  have  developed  a  theory 
of  the  energy  loss  rate  to  acoustic  phonons,  discussed  the  effects  of  static  screening,  and  determined 
a  value  of  the  deformation  potential  constant  by  matching  the  theory  to  our  experimental  data. 
We  have  shown  how  the  temperature  dependence  of  the  energy  loss  rate  can  be  used  to  determine 
the  effectiveness  of  screening.  In  agreement  with  previous  studies,  we  find  that  an  anomalously 
large  deformation  potential  constant  of  16.0  or  U.5eV  for  the  screened  or  unscreened  theory, 
respectively,  is  necessary  to  explain  the  data  ( 1  <>]. 

INTERACTIONS  AND/OR  TECHNOLOGY  TRANSFER; 

Some  of  the  simulation  codes  for  transient  transport  have  been  developed  in  cooperation  with 
the  L.  S.  Army  Electronics  Technology  and  Device  Laboratory.  Ft.  Monmouth. 


PUBLICATIONS  AND  REFERENCES 


JSEP-SPONSORED  PUBLICATIONS: 

[l]  k.  Hess.  T.  k.  Higman.  M.  A.  Emanuel,  and  J.  J.  Coleman.  "New  ultrafast  switching  mechan¬ 
ism  in  semiconductor  heterostructures."  J  Appl.  Phys  .  vol.  60.  no.  10.  pp.  3775-3777.  1986. 
(JSEP  ONR  ARO) 

[2j  J.  Higman  and  k.  Hess.  "Comment  on  the  use  of  the  electron  temperature  concept  for  non¬ 
linear  transport  problems  in  semiconductor  p-n  junctions."  St>lid-State  Electronics .  vol.  29. 
no.  9.  pp.  915-918.  1986.  (JSEP) 

[3]  I.  C.  kizilvalli.  k.  Hess.  G.  J.  Iafrate.  and  D.  Smith.  "Dynamics  of  electron  transfer  between 
two  adjacent  channels  as  calculated  by  an  ensemble  Monte  Carlo  method."  Proc  WAtOS 
Workshop.  Los  Angeles.  CA.  Dec.  1986.  (ARO' JSEP) 


I 

I 

\ 

1 

i 

i 

I 

I 

i 

3 


'I 


a 

& 


[4J  J  P  Leburton.  "Anomalous  current  oscillations  in  semiconductor-insulator-semiconduclor 
structures  and  related  devices."  NATO  Summer  Schml,  Antwerp,  Belgium,  1986.  (JSEP) 

» 

[5]  J.  P  leburton,  "Microscopic  interpretation  of  the  current  oscillations  in  GaAs+  -AlGaAs- 
GaAs~  tunnel  junctions."  Proc.  1 8th  1  C  PS  .  Stockholm.  Sweden.  1986.  (JSEP) 

[6]  S  1  Chuang  and  k  Hess.  "Tunneling-assisted  impact  ionization  for  a  superlattice."  J  Appl 
Pins  vol  61.  no.  4.  Feb.  1987.  (NASA 'JSEP) 

[7]  J  P  I  eburton.  "Space  charge  generation  by  1.0  phonon  ionization  of  neutral  donors  in  GaAs- 
AlGaAs  GaAs  tunnel  junctions."  Bull  Am.  Phys  Soc..  col.  32.  pp.  886.  1987.  (JSEP) 

[8]  T  k  Higman.  J.  M.  Higman.  M.  A.  Emanuel,  k.  Hess,  and  J.  J.  Coleman.  "Theoretical  and  ex¬ 
perimental  analysis  of  the  switching  mechanism  in  heterostructure  hot  electron  diodes."  to  be 
published  in  /  Appl  Phys  ( ARO/ONR/JSEP) 

[9]  1  C  kizilvalli.  k.  Hess,  and  G.  J.  lalrate.  "Electron  transfer  between  adjacent  channels  simu¬ 
lated  by  ensemble  Monte  Carlo  methods."  to  be  published  in  J  Appl.  Phys.  (ARO/JSEP) 

[10]  S  J  Manion.  M.  Artaki.  M.  A.  Emanuel,  J.  J.  Coleman,  and  k.  Hess.  "Electrons-energv-loss 
rates  in  AlGa.As'GaAs  heterostructures  at  low  temperatures."  to  be  published  in  /.  Appl 
Phys  .  (ONR  JSEP) 

PUBLICATIONS  UNDER  OTHER  SPONSORSHIP; 

[ill  k.  Yokoyama  and  k.  Hess.  "Monte  Carlo  study  of  electronic  transport  in 
A 1  ] _xGax As/Ga As  single-well  heterostructures."  Phvs.  Rev.  B.  vol.  33.  no.  8.  pp.  5595-5606. 
1986  (ARO'ONR) 

[12]  T  Wang,  k  Hess,  and  G.  J.  Iafrate.  "Monte  Carlo  simulations  of  hot-electron  spectroscopy  in 
planar-doped  barrier  transistors."  J.  Appl.  Phys..  vol.  59.  no.  6.  pp.  2125-2128.  1986.  (ARO) 

[13]  K.  kim.  k.  kahen.  J.  P.  Leburton.  and  K.  Hess.  "Band-structure  dependence  of  impact  ioniza¬ 
tion  rate  in  GaAs,"  /.  Appl.  Phys.,  vol.  59.  no.  7,  pp.  2595-2596.  1986.  (ONR) 

[14]  G  E.  Stillman.  V.  M.  Robbins,  and  K.  Hess.  "Impact  ionization  in  InP  and  GaAs."  Physica, 
vol.  134B.  pp.  241-246.  1986.  (NSF/ONR) 

[15]  K.  Yokoyama  and  k.  Hess,  "Calculation  of  warm  electron  transport  in  AlGaAs/GaAs  single 
heterostructures  using  a  Monte  Carlo  method."  /.  Appl.  Phvs..  vol.  59.  no.  11.  pp.  3798-3802, 
June  1986.  (NSF) 

[16]  k.  kim  and  K.  Hess.  "Simulations  of  electron  impact  ionization  rate  in  GaAs  in  nonuniform 
electric  fields."  /.  Appl.  Phys.,  vol.  60.  no.  7,  pp.  2626-2629.  Oct.  1986.  (ONR) 

[  1 7]  k.  Hess.  Large  Scale  Computational  Device  Modeling.  Proceedings  of  a  workshop  on  large 
scale  computational  device  modeling,  Apr.  18-19.  1985,  Naperville.  IL.  K.  Hess,  Editor.  1986. 
(NSF) 


-  VaVaVV/.'/aWn'.' 


C.VA  .V.V 


12 


[  1 8 ]  P.  A  Martin.  k.  Hess.  M.  Emanuel,  and  J.  J  Coleman.  "Deep  level  transient  spectroscopy 
studies  of  defects  in  GaAs-AlGaAs  superlattices"  J  Appl.  Ph\s.,  vol.  60.  no.  8.  pp.  2882- 
2885.  Oct.  1986.  (OXR/ARO) 

[  1 9]  I.  C.  Kizilvalli.  k  Hess.  J.  L.  Larson,  and  D.  Widiger,  "Scaling  properties  of  high  electron  mo¬ 
bility  transistors."  IEEE  Trans,  on  Electron  Devices,  vol.  ED-33,  no.  10,  pp.  1427-1433.  Oct. 
1986.  (OXR/ARO) 

[20]  k  B.  kahen  and  J.  P.  Leburton.  "Optical  constants  of  GaAs-AlGaAs  superlattices  and  multi¬ 
ple  quantum  wells."  Phys.  Rev.  B,  vol.  33.  pp.  5465-5472.  1986.  (NASA) 

[2 1  ]  J.  P.  Leburton  and  k.  B.  kahen.  "Modeling  of  the  optical  properties  of  III— V  semiconducting 
compounds  and  related  heterostructures,"  in  Proc.  of  1st  workshop  on  Large  Scale  Computa¬ 
tional  Device  Modeling,  k.  Hess.  Ed..  University  of  Illinois.  Urbana.  IL  61801.  1986,  pp. 
119-130.  (NASA) 

[22]  M.  Wagner  and  J.  P  Leburton.  "A  proposal  for  superstructures  based  high  efficiency  photo- 
voltaics."  Appl.  Phys.  Lett.,  vol.  49,  pp.  886-888,  1986.  (NASA) 

[23]  k.  B.  kahen  and  J.  P.  Leburton.  "Exciton  effects  in  the  index  of  refraction  of  multiple  quan¬ 
tum  wells  and  superlattices."  Appl  Phys.  Lett.,  vol.  49.  pp.  734-736.  1986.  (NASA) 

[24]  D.  Arnold  and  k.  Hess.  "Barrier  height  fluctuations  in  very  small  devices  due  to  the  discrete¬ 
ness  of  the  dopants."  to  be  published  in  J.  Appl.  Phys.  (ONR.  NSF) 

[25]  K.  Hess  and  G.  J.  lafrate.  "Modern  aspects  of  helerojunction  transport  theory."  to  be  pub¬ 
lished. 

[26]  k.  B.  Kahen  and  J.  P.  Leburton.  "Index  of  refraction  of  GaAs-AlGaAs  superlattices  and  mul¬ 
tiple  quantum  wells."  Superlattices  and  Microstructures .  in  press.  (NASA) 


•N-vXk  r.: 


13 


WORK  UNIT  NUMBER  4 

TITLE:  Basic  Studies  of  the  Optical  and  Electronic  Properties  of  Defects  and  Impurities  in  Com¬ 
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SCIENTIFIC  OBJECTIVE: 

The  objective  of  this  research  unit  is  to  contribute  to  our  understanding  of  impurity  incor¬ 
poration  mechanisms,  sources,  and  defects  and  to  improve  our  understanding  of  the  influence  of 
growth  conditions  on  impurities  and  defects  in  semiconductor  materials  that  will  be  important  for 
new  multiple-layer  heterostructure  devices.  It  includes  developing  new  characterization  techniques 
that  will  extend  the  range  of  impurity  concentrations  over  which  quantitative  analysis  is  possible. 

SUMMARY  OF  RESEARCH: 

The  devices  and  multiple  heteroepitaxial  structures  of  interest  for  advanced  electronic  and 
optoelectronic  applications  all  require  the  careful  control  of  alloy  composition  and  doping  levels. 
Some  of  these  devices  require  high-purity  and  low-carrier  concentration  layers  as  well,  as  for 
MODFETs.  HEMTs.  and  quantum  well  avalanche  photodiodes.  We  have  developed  characterization 
techniques  (Fourier  transform  far-infrared  photothermal  ionization  spectroscopy  and  near-infrared 
variable  temperature  photoluminescence)  for  identifying  and  quantifying  the  residual  impurities  in 
high-purity  GaAs.  These  techniques  are  best  suited  to  the  characterization  of  very  high-purity 
materials.  When  the  residual  impurity  concentrations  are  high,  as  is  often  the  case  during  the 
development  of  new  semiconductor  materials,  or  for  small  bandgap  materials,  it  would  be 
extremely  useful  to  have  additional  techniques  which  are  applicable  to  these  cases.  In  this  research, 
another  limitation  of  currently  available  techniques  is  that  it  is  not  presently  possible  to  identify 
donor  species  in  p-type  high-purity  material.  We  proposed  to  develop  high  magnetic  field  (9T). 
resonantly  pumped  high-resolution  Fourier  transform  magneto-photoluminescence  techniques  for 
the  identification  of  residual  donor  and  acceptor  impurities  in  relatively  impure  n-  and  p-type 
GaAs.  These  techniques  will  also  be  applied  to  InP.  AlGaAs.  InGaAs.  and  other  related  alloys. 
This  work  will  be  critical  for  the  development  of  high-puritv  MBE  and  MOCVD  source  materials 
for  new  compound  semiconductor  materials. 
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The  Fourier  Transform  photoluminescence  technique,  with  the  use  of  very  high  magnetic 
fields,  will  also  permit  photoluminescence  measurements  to  he  extended  to  small  bandgap  materials 
such  as  InGaAs.  HgCdTe.  and  related  semiconductors.  These  measurements  are  not  very  effective 
using  ordinary  dispersive  techniques. 

During  this  reporting  period  we  have  achieved  identification  of  residual  donors  in  high-puritv 
epitaxial  GaAs  by  magneto-photoluminescence  (MPI.)  measurements  at  a  high  magnetic  field 
(9.0  T)  with  selective  excitation.  The  high-purity  epitaxial  layers  used  for  this  identification  were 
grown  by  various  techniques  which  have  characteristic  residual  impurities.  The  results  of  77  k 
Hall  effect  measurements  on  these  samples  confirm  the  high  purity  of  t he  layers.  The  different 
donor  species  have  been  detected  in  these  layers  from  well-resohed  "tw o -electron"  satellites  in  the 
MPI  spectra.  Photothermal  ionization  spectroscopic  il'TISi  measurements  were  also  made  on  the 
same  layers  to  identity  the  donor  species  present.  The  MPI  anil  l’TIS  measurements  have  been 
correlated  to  identify  the  "two-electron"  satellite  peaks  in  the  MPI  spectra  corresponding  to  the 
different  donor  species. 

Magnetophotoluminescence  measurements  were  made  on  the  samples  mounted  strain  free  in 
pumped  liquid  helium,  at  a  high  magnetic  field  (9.0  T)  in  a  superconducting  solenoid  in  the  Fara¬ 
day  configuration.  An  infrared  tunable  dye  laser  with  a  power  of  (20-30  mW)  was  used  to  pump 
the  luminescence.  Styryl  9  M  dye  was  used  to  cover  the  range  of  excitation  energies  slightly  below 
and  above  the  bandgap  of  GaAs.  The  dye  laser  was  optically  pumped  by  -  2-4  watts  of  all  lines 
Irom  an  argon  ion  laser  An  arrangement  of  lenses  was  used  to  delocus  the  laser  beam  onto  the 
sample.  The  luminescence  was  collected  and  focussed  onto  a  four-slit  1.0  m  double  spectrometer 
and  detected  by  a  ihermoelectrically  cooled  GaAs  photomultiplier  tube  using  the  photon  counting 
technique. 

Figure  1  shows  the  typical  PL  spectrum  ol  exciton  recombination  in  high-puritv  epitaxial 
GaAs  (sample  E)  grown  by  AsCl,— \  PE  (capor  phase  epitaxy  ).  The  spectrum  was  recorded  at 
1.7  k.  low  excitation  intensity  and  a  high  magnetic  field  (9  0  T).  The  dye  laser  output  was  tuned 
close  to  the  bandgap  of  GaAs.  The  highest  energy  peak  labeled  FE'  is  the  free  exciton  recombina¬ 
tion.  The  large  peaks  in  the  high-energv  region  between  1.5 ISO  eV  and  1.5193  eV  are  the  princi¬ 
pal  lines  of  the  donor-bound  exciton  (D°.\)  emission.  The  principal  lines  labeled  (D°.X)ls  are  tran¬ 
sitions  in  which  the  initial  state  is  the  ground  state  or  an  excited  slate  of  the  (D°.X)  complex,  and 
the  final  state  is  the  ground  state  of  the  donor.  The  peaks  in  the  low-energy  region  between 
1.510  eV'  and  1.515  eV  are  the  "two-electron"  transitions.  These  transitions  differ  in  energy  from 
the  principal  lines  by  the  energy  (above  the  donor  ground  stale)  of  the  final  state  of  the  transition 
in  which  the  donor  is  left  in  its  first  excited  state.  At  high  magnetic  fields,  "two-electron"  transi¬ 
tions  which  leave  the  donors  in  the  2p„  and  2p_j  terminal  excited  states  are  observed  to  be  the 
strongest.  Moreover,  only  two  transitions  to  the  2p_|  magnetic  substate  and  one  transition  to  the 
2p„  magnetic  substate  are  observed  to  have  significant  intensities  in  the  “two-electron"  satellite  spec¬ 
trum.  The  "two-electron"  transition  is  not  observed  to  be  strong  for  every  single  principal  line  due 
to  selection  rules. 

The  high-resolution  spectra  of  "two-electron"  transitions  were  obtained  by  resonant  excitation 
of  one  of  the  strong  principal  lines  of  (D°,Xj.  Figure  2  shows  the  "two-electron"  satellite  spectra 
for  the  samples  E  and  A.  The  top  spectrum  in  Figure  2  lor  this  sample  was  recorded  with  the  dve 
laser  output  tuned  to  the  line  b'  in  Figure  F  The  bottom  spectrum  was  recorded  for  sample  A 
grown  by  MOCVD  ( metalorganic  chemical  vapor  deposition)  with  the  dve  laser  output  tuned  into 
resonance  with  one  of  the  strong  principal  lines,  e.g..  the  line  V  in  this  case.  To  identify  the  "two- 
electron"  MPL  peaks  corresponding  to  particular  donor  species  the  results  of  MPL  measurements  at 
different  magnetic  fields  have  been  correlated  with  the  results  of  PT1S  measurements  on  the  same 
samples.  This  correlation  along  with  the  results  of  resonance  excitation  has  been  used  to  relate  the 
"two-electron"  satellite  peaks  to  the  principal  lines.  The  principal  line  labeled  b'  in  Figure  1  and 
the  low-energy  strong  line  labeled  '3  in  the  "two-electron"  transition  to  the  2p_,  donor  magnetic 
substate  are  found  to  be  related  to  each  other  in  that  the  same  initial  state  is  invoked  in  both  tran¬ 
sitions.  The  correlation  suggests  that  the  single  strong  line  in  the  "two-electron"  transition  to  the 


■  «  ^*«  ■**«  f  W* 


.V-VjLSJC. 


>A. 


"Two-electron"  spectra  for  two  GaAssamples  recorded  at  a  magnetic  field  H  »  9.0  T  with 
the  resonant  excitation  of  one  of  the  strong  principal  lines.  Si  and  Ge  donors  are  identified 
in  the  samples  E  and  A.  respectively. 
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2p,,  donor  magnetic  substate  is  a  "two-electron"  replica  of  the  principal  line  'b‘  or  a  line  close  to  b 
The  energy  separation  between  the  principal  line  b'  and  the  corresponding  "two-electron"  satellite 
peak  (2p_!  case)  labeled  /3  gives  the  Is — 2p_ ,  energy  for  the  donors.  The  principal  line  correspond¬ 
ing  to  the  high-energy  strong  line  in  the  "two-electron"  transition  to  the  2p_,  donor  magnetic  sub¬ 
state  is  not  observed.  This  unobserved  principal  line,  which  is  very  close  to  the  line 'd  ,  must  be  an 
unfavored  transition  because  of  the  influence  of  the  selection  rules.  The  Is— 2 p_ i  energies  can  be 
used  to  identify  donors  from  the  MPI.  spectra.  The  donors  identified  in  samples  E  and  A  from 
MPI  spectra  are  Si  and  Ge.  respectively  The  dominant  Ge  donor  in  sample  A  has  been  identified 
from  the  correlation  ol  PT1S  measurements  with  MPI  measurements.  The  identification  of  Si 
donors  in  the  sample  F  is  also  consistent  with  the  correlation.  The  Is— 2p_]  energies  for  Si  and  Ge 
donors  measured  from  the  MPT  spectra  tor  samples  1  and  A  are  4.64  meV  and  4.87  meV.  respec¬ 
tively.  at  a  magnetic  field  of  9.0  T 

Figure  3  shows  the  "two-electron"  satellite  spectra  for  samples  B.  C.  and  D  recorded  at  a  mag¬ 
netic  field  of  9.0  T  with  the  dye  laser  tuned  to  a  strong  principal  line  for  samples  C  and  D.  The 
spectrum  for  sample  B  was  recorded  with  above  bandgap  excitation.  The  samples  are  very  high- 
purity  GaAs  epitaxial  layers  grown  by  liquid  phase  epitaxy  (LPE)  and  molecular  beam  epitaxy 
(MBE).  S  is  identified  as  a  dominant  donor  in  the  LPE  GaAs  sample  B.  The  donors  identified  in  the 
MBE  GaAs  sample  C  (lightly  Si-doped)  are  S  and  Si.  Also.  S.  Sn  (Se).  and  Te  donors  are  identified 
in  the  third  LPE  GaAs  sample  D  The  Is— 2p_[  energies  for  Si  and  S  donors  measured  from  MPL 
spectra  are  4.64  meV  and  4.73  meV.  respectively,  at  a  magnetic  field  of  9.0  T. 

The  energy  separation  of  the  peaks  due  to  different  donors  in  the  MPL  spectra  of  the  "two- 
electron"  satellite  gives  the  central  cell  shifts  between  the  Is  (ground)  states  of  the  donors.  The 
central  cell  shift  between  S  and  Si  measured  from  the  "two-electron"  satellite  spectrum  for  the 
MBE  GaAs  sample  is  0.12  meV  at  a  magnetic  field  of  9.0  T  The  central  cell  shifts  between  S  and 
Sn  (or  Se)  donors  and  between  S  and  Te  donors  measured  from  the  satellite  spectrum  for  the  sam¬ 
ple  D  are  0  09  meV  and  0.18  meV.  respectively,  at  a  magnetic  field  of  9.0  T.  The  central  cell 
shifts  between  donors,  measured  from  the  satellite  spectra  for  samples  in  which  more  than  one 
donor  species  have  been  detected,  are  little  larger  than  what  should  be  expected  based  on  the  corre¬ 
lation  of  PTIS  measurements  with  MPL  measurements. 

All  the  major  donors.  S.  Si.  Ge.  Sn  (or  Se).  and  Te  are  identified  by  magneto¬ 
photoluminescence  measurements  on  the  GaAs  epitaxial  layers  grown  by  various  techniques.  These 
results  manifest  good  correlation  with  the  photothermal  ionization  spectroscopic  measurements  on 
the  same  epitaxial  layers. 

During  the  second  year  of  this  work,  photoluminescence  measurements  by  Fourier  Transform 
spectroscopy  will  be  developed.  This  should  result  in  faster,  higher  resolution  measurements  and 
when  developed  will  permit  the  photoluminescence  study  of  narrow  bandgap  materials  using  the 
same  technique 
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"Two-electron"  spectra  for  three  GaAs  samples  recorded  at  a  magnetic  field  of  9.0  T.  The 
principal  lines  occur  at  nearly  the  same  energies  for  these  samples.  Spectra  for  LPE  and 
MBE  GaAs  samples  show  that  the  peaks  due  to  S  donors  match  in  energy  position  and  the 
peaks  corresponding  to  Sn  (or  Se)  and  Te  donors  lie  on  both  sides  of  the  Si-peak. 
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SCIENTIFIC  OBJECTIVE: 

The  objective  of  this  program  is  to  extend  to  electronic  devices  the  enormous  impact  that 
metalorgamc  chemical  vapor  deposition  (MOCVD).  as  a  sophisticated  epitaxial  growth  method,  has 
had  on  optical  device  research.  This  involves  fundamental  studies  of  the  MOCVD  process  itself  for 
electronic  materials,  studies  of  the  electronic  properties  of  heterostructure  electronic  materials,  and 
studies  of  devices  made  from  these  materials.  Three  specific  areas  of  interest  for  this  research  are 
( 1 )  development  of  new  reaction  chamber  designs  and  growth  regimes  for  MOCVD-grown  heteros¬ 
tructure  electronic  materials.  (2)  continuation  of  structural  and  materials  analysis  including  deep- 
level  transient  spectroscopy  (DLTS)  and  Shubnikov-de  Haas  measurements  of  superlattice  struc¬ 
tures.  and  (3)  development  of  MOCVD-grown  real  space  transferred  electron  devices  and  other 
electronic  devices. 
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SUMMARY  OF  RESEARCH: 

In  the  past  year,  we  have  developed  the  heterostructure  hot  electron  diode  (H’ED).  a  new 
two-terminal  device  that  exhibits  S-shaped  negative  differential  resistance  (XDR)  due  to  a  field- 
dependent  transition  between  the  current  conduction  modes  of  tunneling  and  thermionic  emission 
of  hot  electrons  in  a  two-layer  AlGaAs-GaAs  heterostructure.  Experimental  results  have  been 
obtained  for  various  single  period  and  multiple  period  H3ED  structures  fabricated  from  wafers 
grown  by  metalorgamc  chemical  vapor  deposition  (MOCVD).  We  have  shown,  using  an  analvtic 
theory,  that  the  field-switching  mechanism  previously  proposed  is  in  fact  responsible  for  the  NDR 
observed  in  the  heterostructure  hot  electron  diode.  In  addition  to  the  tunneling  and  thermionic 
emission  of  hot  electrons,  the  accumulation  of  electrons  at  the  heterointerface  is  identified  as  an 
important  mechanism  in  the  operation  of  the  device.  Experimentally,  we  verify  that  the  basic  phy¬ 
sical  requirements  for  switching  predicted  by  the  theory  are  correct  and  further  we  show  that  the 
potentially  detrimental  effect  of  impact  ionization  can  be  minimized  by  proper  design.  Preliminary 
investigations  of  the  H2ED  on  nonoptimized  structures  have  resulted  in  test-fixture-limited  oscilla¬ 
tion  at  frequencies  greater  than  17  GHz.  These  results  indicate  that,  as  reported  previously,  the 
H'ED  should  be  capable  of  transit  time-limited,  high-speed  performance. 
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SCIENTIFIC  OBJECTIVE: 

Objectives  are  twofold.  The  first  is  to  characterize  and  understand  the  optical  and  vibrational 
properties  of  novel  semiconductor  superlattices  and  multiple  quantum  wells.  The  second  is  to 
studs  and  understand  the  dynamics  of  photopumped  carriers  in  quantum  wells. 


SUMMARY  OF  RESEARCH: 


fV  The  JSFP-supported  portion  of  our  work  on  Raman  and  other  optical  properties  of  semicon- 

H  ductor  quantum  wells  and  superlattices  has  consisted  of  two  parts:  (1)  a  Raman  study  of  pho- 

>k  nons  in  a  GaAs-AlAs  superlattice  as  the  two  components  diffuse  into  each  other  upon  high  tempera- 

v  ture  annealing:  and  (2)  setting  up  equipment  for  time-resolved  absorption  measurements  of  photo- 

pumped  GaAs  quantum  wells. 


Properties  of  Phonons  in  a  GaAs-AlAs  Superlattice  upon  Annealing 

Since  joining  the  JSEP  program,  we  have  used  Raman  scattering  to  study  the  effects  of  super- 
lattice  layering  on  the  acoustical  and  optical  phonons  in  semiconductors.  Our  earlier  work  has 
shown  conclusively  that  very  simple  models  apply  to  phonons  having  wavevector  q  normal  to  the 
interfaces.  The  longitudinal  acoustical  (LA)  modes  propagate  freely  through  the  interfaces, 
suffering  only  relatively  weak  Bragg  reflection,  and  have  dispersion  curves  given  by  minizone  fold¬ 
ing  of  the  virtual-crystal  dispersion  curve.  The  longitudinal  optic  (LO)  modes  are  strongly 
confined  into  either  Ga.As-rich  or  AlAs-rich  regions  and  show  little  dispersion  with  q. .  We  have 
now  completed  an  extensive  study  of  how  these  modes  change  when  the  concentration  profile  of  the 
superlattice  is  changed  by  high-temperature  annealing.  That  profile  was  measured  using  x-ray 
diffraction.  The  folded  acoustic  phonons  did  not  shift  in  frequency  but  did  change  intensity  upon 
annealing.  The  intensity  of  the  two  lowest  doublets  followed  that  of  the  two  smallest  superlattice 
diffraction  peaks,  in  agreement  with  our  photoelastic  model  for  coupling  of  light  to  acoustic  pho¬ 
nons  in  a  superlattice.  The  third-order  doublet  was  anomalously  strong  due.  we  believe,  to  a 
resonant  enhancement  of  the  photoelastic  coefficient  in  Al-rich  portions  of  the  superlattice.  We  fol¬ 
lowed  up  to  the  ninth  order  of  confined  LO  phonon  modes.  With  annealing,  they  shift  down  in  fre¬ 
quency  due  to  an  effective  narrowing  of  the  GaAs  phonon  quantum  wells  due  to  diffusion  of  Al. 
This  can  be  quantitatively  explained  using  an  effective  mass  model  for  the  phonons.  The  Al  con¬ 
centration  profile  generates  the  "potential  energy"  term  in  the  wave  equation  through  the  known 
dependence  of  LO  frequency  on  concentration.  Known  phonon  dispersion  for  GaAs  gives  the 
"effective  mass"  term.  These  results  show  that  a  simple  model  exists  for  LO  phonons  in  quantum 
wells  and  superlattices  having  a  general  concentration  profile. 
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Instrumentation  for  Time-Resolved  Absorption  in  Quantum  Wells 

In  the  Laser  Laboratory  of  the  Materials  Research  Laboratory,  we  have  set  up  equipment  that 
allows  us  to  measure  the  change  in  absorption  in  quantum  well  samples  produced  by  pumping  car¬ 
riers  into  the  conduction  subbands.  With  JSFP  and  other  support,  we  have  ordered  a  femtosecond, 
cav  it y '-dumped  dye  laser  that  will  extend  our  capability  for  time-resolved  absorption  measure¬ 
ments  into  the  sub-picosecond  regime. 

Molecular  Beam  Epitaxy 

During  this  period  we  have  concentrated  our  efforts  on  understanding  the  material  quality  of 
1  \!.(  »a  'As  ( iu As  and  (In.Ga)As  GaAs  multi-quantum  w  ell  ( MQW)  structures  on  compound  and 
elemental  substrates  grown  bv  molecular  beam  epitaxy  (MBL).  Room  temperature  photoreflectance 
measurements  were  used  in  these  studies.  From  the  study  ol  ( Al.GaJAs  GaAs  MQWs  [9]  grown  on 
n±  GaAs  substrates,  we  have  observed  transitions  involving  the  so  called  "unconfined"  states, 
which  have  received  less  attention  so  far.  The  study  of  a  low-barrier  multiple  quantum-well 
structure  has  allowed  us  to  observe  many  transitions  as  far  as  200  meV  beyond  the  barrier  gap. 

The  other  photoreflectance  studies  [4.10]  on  the  ( Al.Ga )As/GaAs  and  GaAs/(In,Ga)As  MQWs 
grown  on  Si  and  Ge  have  demonstrated  that  a  variety  of  III- V  compound  devices  can  be  grown  on 
elemental  substrates. 

In  a  separate  investigation  we  have  studied  the  effect  of  annealing  on  the  dislocation  reduction 
near  the  interface  of  the  GaAs/Si  and  in  the  bulk  of  Ga.As  layers  [3],  We  have  demonstrated  that 
ex-situ  and  in-situ  annealing  at  higher  temperatures  reduces  the  threading  dislocations  consider¬ 
ably.  This  will  aid  in  growing  high  quality  active  layers  for  high  speed  applications  on  Si  sub¬ 
strates. 

Finally,  we  have  studied  the  mobility  degradation  in  GaAs  and  (In.Ga)As  [5],  GaP  [6]  and 
( Al.GaLAs  GaAs  [7]  structures  at  high  electric  fields.  We  have  shown  that  the  asymmetry  in  the 
distribution  function  of  the  electrons  is  the  cause  of  the  mobility  degradation.  Furthermore,  the 
mobility  of  high-mobility  materials  degrades  faster  with  the  electric  fields  than  low-mobility 
materials. 
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SCIENTIFIC  OBJECTIVE: 

The  continuing  evolution  of  integrated  circuit  processing  technology  toward  ultra-small 
feature  sizes  and  new  compound  semiconductor  materials  and  devices  demand  new  computer-aided 
design  (CAD)  tools  and  design  methods.  One  of  the  basic  issues  is  to  develop  efficient  circuit 
models  for  new  devices  such  as  the  High-Electron  Mobility  Transistor  (HEMT)  and  important 
parasitics  of  interconnects  in  high-speed  integrated  circuits.  With  such  models,  integrated  circuits 
can  be  simulated  accurately.  Apart  from  modeling,  problems  related  to  the  chip  yield  and  reliabil¬ 
ity.  design  verification  of  L'ltra-/Very-Large  Scale  Integrated  (L’/VLSI)  circuits  by  using  novel 
parallel  computation  algorithms  and  hardwares,  and  the  chip-performance  optimization  in  view  of 
yield  and  reliability  are  becoming  issues  of  increasing  importance. 

The  objectives  of  our  research  in  this  unit  are: 

( 1 )  to  develop  a  computationally  efficient  model  for  the  HEMT  and  interconnection  parasitics  so 
that  high-speed  integrated  circuits  can  be  simulated  accurately  in  short  time: 

(2)  to  develop  rigorous  design  optimization  techniques  for  yield  and  reliability  enhancement:  and 

(3)  to  develop  accelerated  circuit  simulation  methods  for  L7VLSI  circuits  by  exploiting  hierar¬ 
chy  and  parallelism. 


SUMMARY  OF  RESEARCH: 

This  past  year,  a  new  physical  model  has  been  developed  and  implemented  in  the  Illinois 
Simulator  for  Modeling  of  Integrated-circuit  Level  Elements  (iSMILE)  program,  which  was 
developed  under  the  JSEP  sponsorship.  This  model  is  charge-based  and  uses  a  novel  analytical 
approximation  of  channel-charge  versus  gate-voltage  characteristics,  which  renders  a  closed-form 
solution  for  the  channel  current.  The  model  equations  for  the  current  and  charge  have  continuous 
derivatives  and  hence  are  well-suited  for  circuit  simulation.  With  accurate  dc  characterization  of 
both  enhancement  and  depletion  devices,  we  have  been  able  to  predict  ring  oscillator  circuit  delays 
with  less  than  10%  error.  Device  samples  have  been  obtained  from  both  AT&T  Bell  Laboratories 


(Dr.  S.  S  Pei)  and  the  Coordinated  Science  Laboratory  (Prof.  H.  Morkoc.).  The  effects  of  parasitic 
capacitances  and  resistances  have  been  investigated  extensively  to  determine  optimal  device 
geometries  for  maximum  circuit  speed.  By  using  our  analytical  equations  for  the  gate  fringing 
capacitance  and  the  interelectrode  capacitances,  we  can  predict  the  optimal  source-to-gate  spacing 
w  hich  seems  to  match  well  with  the  current  state-of-the-  art  devices. 

This  research  will  be  continued  in  the  coming  year  to  extend  the  domain  of  the  model  to 
small-signal  analysis  of  high-power  microwave  circuits.  This  would  include  the  improvement  and 
v  erification  of  the  model  with  S-parameter  measurements  at  diff  erent  dc  bias  points  and  at  many 
ditlerent  frequencies.  Also  much  more  w'ork  needs  to  be  done  on  the  modeling  of  back-gate  biasing 
effects  and  the  gate-diode  capacitance.  Statistical  analysis  is  also  required  to  determine  the  effects 
ot  various  process  deviations  on  the  device  yield. 

for  the  improvement  of  yield  and  reliability,  we  have  developed  a  new  hierarchical  timing 
optimizer  iJADE  for  MOS  VLSI  circuits.  iJADE  uses  a  switch-  level  timing  simulator  for  accurate 
delay  calculation  and  a  rule-based  expert  system  to  optimally  determine  individual  transistor  sizes 
for  speed-area  trade-offs.  iJADE  takes  in  the  SPICE  input  descriptions  for  the  circuit  to  be  optim¬ 
ized.  partitions  the  circuit  hierarchically  to  save  the  memory  space,  detects  critical  paths,  and  then 
finds  optimal  transistor  sizes  to  meet  the  timing  specifications.  iJADE  is  written  in  Franz  LISP 
language  and  at  present  consists  of  about  10.000  lines  of  code.  It  runs  on  the  VAX-1 1/780  under 
the  L'NIX  operating  system.  Several  circuits  have  been  tested  with  iJADE.  and  it  has  been  demon¬ 
strated  that  iJADE  could  find  a  better  circuit  construct  than  a  design-in-mask  from  industry  con¬ 
taining  more  than  140  transistors.  Also.  iJADE  was  able  to  eliminate  all  the  glitches  in  the  industry 
design. 

To  make  our  circuit  optimization  system  more  versatile,  this  work  needs  to  be  continued  in 
the  coming  year  in  order  to  test  larger  and  more  difficult  examples.  It  is  desirable  to  include  a 
self-learning  capability  so  tlie  rules  for  circuit  optimization  can  evolve  automatically  as  iJADE  is 
exercised  for  more  difficult  circuits. 

Multiprocessors  offer  a  way  to  reduce  the  time  required  to  perform  circuit  simulations  of 
large  circuits.  To  obtain  this  speed  advantage,  the  simulation  algorithm  must  have  sufficient  con¬ 
currency  to  keep  the  processors  usefully  employed  without  introducing  excessive  memory  conten¬ 
tion  and  scheduling  overhead.  Future  multiprocessors  with  larger  numbers  of  processors,  such  as 
the  Cedar  machine  which  is  to  consist  of  multiple  8-processor  Alliant  FX/8  computers,  will  require 
algorithms  with  higher  degrees  of  concurrency  in  order  to  achieve  further  speed  improvements. 

We  have  studied  alternative  strategies  for  employing  waveform  relaxation  (WR)  algorithms 
to  perform  circuit  simulations  of  digital  MOS  circuits  on  multiprocessors.  WR  has  a  potential 
speed  advantage  due  to  its  multirate  integration  properties,  and  it  exhibits  natural  concurrency  due 
to  the  inherent  decoupling  of  the  circuit  into  independently  solved  subcircuits. 

To  study  the  performance  of  alternative  WR  algorithms  on  an  actual  multiprocessor,  we 
developed  a  WR  circuit  simulator  to  run  on  an  8-processor  Alliant  FX/8  computer.  Both  the 
Gauss-Seidel  and  Gauss-Jacobi  methods  are  implemented.  In  the  current  version  of  the  program, 
tasks  are  defined  only  at  the  time-window  level.  Circuits  with  up  to  805  transistors  have  been 
simulated.  The  results  demonstrate  that  when  the  number  of  processors  is  sufficiently  large,  the 
Gauss-Jacobi  method  generally  outperforms  the  Gauss-Seidel  method.  The  opposite  is  the  case 
when  the  number  of  processors  is  small.  The  break-even  point  between  the  two  methods  is  a  func¬ 
tion  of  the  properties  of  the  circuit  being  simulated. 

To  determine  whether  to  use  the  Gauss-Seidel  or  Gauss-Jacobi  method  for  a  given  circuit  on  a 
given  number  of  processors,  we  have  developed  a  computationally  inexpensive  algorithm  which 
estimates  the  break-even  point  between  the  two  methods.  The  break-even  point  calculation  uses  an 
estimate  of  concurrency  based  on  the  task  precedence  graph  and  task  CPU  times  which  are  assumed 
to  be  proportional  to  subcircuit  size. 


The  time-point  pipelining  strategy  was  first  reported  in  the  literature  as  a  way  to  retain  the 
fast  convergence  properties  of  the  Gauss-Seidel  method  while  at  the  same  lime  improving  its  con¬ 
currency  by  synchronizing  tasks  at  the  time-point  level  rather  than  the  time-window  level.  We 
have  generalized  the  time-point  pipelining  algorithm  and  applied  it  to  the  (iauss-Jacobi  method  In 
this  algorithm,  all  subcircuits  in  all  iterations  can  be  solved  concurrently,  with  only  a  slight  delay 
from  the  start  of  one  iteration  to  the  start  of  the  next  iteration.  The  high  degree  of  parallelism  of 
this  method  is  of  particular  interest  as  multiprocessors  with  larger  numbers  of  processors  become 
available. 

A  parallel  simulation  time  estimation  program — FAR ASITE — was  developed  to  estimate  the 
run  limes  of  WR  simulations  performed  on  any  specified  number  ol  processors.  PARASITE  uses 
('PL  times  of  individual  tasks  measured  during  an  actual  simulation  run  on  a  uniprocessor  It  then 
simulates  the  assignment  of  tasks  to  processors  and  estimates  the  overall  multiprocessor  C'PL  lime 
neglecting  multiprocessing  overhead  The  resuits  show  that  the  Gauss-Jacobi  method  with  lime- 
point  pipelining  can  outperform  the  Gauss-Seidel  method  with  time-point  pipelining  if  the  number 
of  processors  is  sufficiently  large  and  if  multiprocessing  overhead  can  be  kept  sufficiently  small. 

We  are  currently  implementing  the  generalized  time-point  pipelining  algorithm  in  our  mul¬ 
tiprocessor  WR  program.  The  relationship  between  window  size  and  time-point  pipelining  perfor¬ 
mance  will  be  investigated.  We  also  plan  to  study  the  use  of  parallel  model  evaluation  in  conjunc¬ 
tion  with  the  other  strategies  to  alleviate  bottlenecks  which  arise  due  to  large  subcircuits.  The  pro¬ 
gram  will  be  adapted  to  run  on  the  larger  Cedar  multiprocessor  when  it  becomes  available. 
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SCIENTIFIC  OBJECTIVE: 

The  goal  of  this  research  is  to  understand  the  basic  mechanism  of  charge  density  wave  (CDW) 
transport  in  the  quasi-one-dimensiona!  metals  NbSea  .  TaSj  .  and  similar  materials.  The  motion  of 
CDWs  represents  the  only  known  example  of  current  flow  bv  a  moving  quantum  ground  state, 
apart  from  superconductivity  John  Bardeen  has  proposed  that  CDW  acceleration  in  an  electric 
field  takes  place  via  a  Zener  tunneling  process  [9]  and  that  the  ac  response  can  be  predicted  by 
applying  the  photon-assisted  tunneling  theory  originally  developed  for  superconductor-insulator- 
superconductor  (SIS)  millimeter  wave  mixers  (to],  If  this  interpretation  is  correct.  CDW  materials 
could  become  important  systems  for  use  as  detectors  and  mixers,  as  well  as  for  basic  studies  of 
macroscopic  quantum  tunneling  An  extensive  experimental  program  is  conducted  in  our  labora¬ 
tory  to  characterize  the  dynamics  of  CDW  motion  and  to  test  the  predictions  of  the  tunneling 
theory. 

SUMMARY  OF  RESEARCH; 

The  past  year  has  been  unusually  productive,  with  major  progress  in  the  following 
areas:  (1)  understanding  the  nature  of  the  low-frequency  dielectric  relaxation  behavior  observed 
in  CDW  materials  below  the  dc  threshold  field.  (2)  characterizing  the  effective  pinning  potential 
that  gives  rise  to  current  oscillations  and  ac-dc  interference  phenomena  for  CDW  motion  above  the 
dc  threshold,  and  (3)  the  beginnings  of  a  microscopic  theory  of  CDW  dynamics  that  demonstrates 
why  tunneling  is  required  for  CDW  acceleration. 

Our  work  on  dielectric  relaxation  [4]  has  resolved  a  major  controversy  and  led  to  new  under¬ 
standing.  Experiments  at  Bell  Laboratories  have  characterized  a  dielectric  relaxation  frequency  in 
KdjMoOj.  (TaSe^I.  and  TaS3  which  decreases  with  an  Arrhenius  behavior.  a>0  cc  e  T°t.  at  low  tem¬ 
peratures.  Workers  at  UCLA  saw  no  evidence  for  this  phenomenon  in  their  low-frequency  ac  con¬ 
ductivity  measurements.  Our  experiments  confirmed  the  existence  of  this  effect  and  showed  that 
the  temperature-activated  behavior  was  the  result  of  normal  electron  screening  of  CDW  charge 
fluctuations  associated  with  the  polarization. 

The  work  on  current  oscillations  and  ac-dc  interference  which  we  performed  during  the  past 
year  [l, 2.5. 6. 7]  has  produced  some  spectacular  results.  By  careful  growth  and  mounting  of  \bSe3 
crystals,  we  have  been  able  to  demonstrate  extraordinary  coherence  for  CDW  motion  throughout 
the  entire  sample  volume— at  least  an  order  of  magnitude  greater  coherence  than  previous  work. 
The  results  of  our  experiments  demonstrate  that  the  CDW  experiences  a  highly  nonsinusoidal 
effective  potential  when  moving  above  threshold,  and  that  the  magnitude  of  this  pinning  potential 
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remains  constant  up  to  the  highest  fields  F  —  1(X)  ET  that  can  be  experimentally  applied.  These 
results  are  most  important  because  they  conclusively  invalidate  the  classical  deformable  model  of 
m\V  d\  namics  that  has  served  as  the  paradigm  for  most  work  in  this  field  over  the  past  five  years 

[5]  This  lends  great  support  to  the  tunneling  interpretation,  am)  it  is  not  now  clear  whether  any 
classical  model  can  be  brought  into  line  w  ith  all  our  experimental  results. 

On  the  theory  front,  we  have  written  a  very  long  paper  [N]  which  derives  John  Bardeen's  ori¬ 
ginal  tunneling  model  [9]  by  employing  some  rather  exotic  methods  in  one-dimensional  quantum 
field  theorv  This  work  clarifies  the  basic  arguments  for  C1)\V  tunneling  and  show's  how  phase 
coherence  car.  be  realized  bv  analogy  to  the  Josephson  effect. 
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SCIENTIFIC  OBJECTIVE: 

The  overall  objective  of  this  research  is  to  investigate  fundamental  plasma  and  chemistry 
processes  in  a  type  of  magnetically  confined  discharge  that  has  promise  for  reducing  radiation  dam¬ 
age  during  device  processing  by  dry  etching. 

Program  Organization 

This  program  is  organized  around  three  specific  objectives. 

1  )  To  investigate  the  basic  plasma  behavior  in  rf  driven  cylindrical  magnetrons  and  the  resultant 
plasma-surface  reactions. 

2 )  To  explore  the  plasma  chemistry  that  is  operative  in  low  pressure  rf-driven  magnetron  etch¬ 
ing  plasmas,  beginning  with  an  examination  of  CH.,— H;  and  CF4— 03  plasmas. 

3)  To  investigate  magnetron  reactive  ion  etching  of  Si  with  particular  emphasis  on  radiation 
damage  and  on  the  basic  surface  chemistry  that  establishes  selectivity  and  anisotropic  etching. 


SUMMARY  OF  RESEARCH: 

The  research  is  proceeding  as  described  in  our  original  proposal.  There  have  been  no  changes 
in  the  overall  technical  objectives  (described  above),  the  direction,  the  participants  or  the  resource 
allocation. 

A  plasma  discharge  apparatus  has  been  designed,  built  and  placed  in  operation,  which  is  capa¬ 
ble  of  being  operated  with  cylindrical  magnetion  electrodes  of  the  type  used  for  reactive  ion  etch¬ 
ing.  The  apparatus  can  accommodate  magnetron  electrodes  of  various  diameters  and  is  instru¬ 
mented  for  making:  ( 1 )  electrostatic  probe  measurements  of  the  plasma  and  floating  potentials  and 
the  electron  density  and  energy.  (2)  retarding  grid  electrostatic  analyzer  measurements  of  energies 
of  the  ions  bombarding  the  electrode  and  wall  surfaces.  (3)  optical  emission  spectroscopy  measure¬ 
ments  of  the  electron  excitation  behavior  in  the  plasma,  and  (4)  mass  spectroscopic  measurements 
of  the  species  produced  in  the  plasma. 

Fundamental  studies  of  the  basic  plasma  behavior  in  rf-driven  cylindrical  magnetron 
discharges  have  been  completed.  Argon  and  He  were  used  as  working  gases.  In  addition  to  provid¬ 
ing  fundamental  information,  which  is  being  prepared  for  publication  [l-6],  these  studies  have  pro¬ 
vided  guidelines  for  the  reactive  ion  etching  studies  that  are  presently  underway.  In  particular,  the 
investigation  has  shown  that  the  differences  in  current-voltage  behavior  between  dc  and  rf  driven 
magnetrons,  and  the  differences  in  the  tendency  to  develop  a  dc  electrode  bias  that  is  seen  when 
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comparing  rf  driven  cylindrical  magnetrons  with  rf  driven  planar  diode  discharges  of  the  type 
commonly  used  in  plasma  etching  can  be  explained  in  terms  of  the  electron  transport  across  the 
magnetic  field  The  dc  bias  behavior  is  significant  in  the  context  of  plasma  etching  because  it  affects 
the  energy  distribution  of  ions  which  are  incident  on  the  various  surfaces  in  contact  with  the 
plasma  Because  the  dc  bias  is  dependent  on  the  electron  transport  across  the  magnetic  field,  it  can 
be  influenced  by  both  the  magnetic  field  strength  and  its  configuration.  The  ion  energy  distribu¬ 
tions  that  were  measured  with  the  retarding  grid  electrostatic  analyzers  in  our  experiments,  and 
which  contributed  to  the  above  conclusions,  were  consistent  with  the  measured  plasma  potential 
distribution  one.  w  ith  t lie  sputtering  which  occurred  w  hen  Ar  was  used  as  a  working  gas. 

Theoretical  computer  studies  have  been  completed  to  assist  in  the  data  interpretation.  These 
studies.  ,v  in.. h  apple  tv'  both  magnetrons  and  plasmas  diodes,  are  being  prepared  for  publication  [6 ]. 
They  show  the  importance  ot  the  discharge  frequency  and  explain  much  of  the  data  that  has  been 
reported  in  the  literature.  Similar  studies  relating  to  non-magnetron  discharges  are  underway  at 
the  University  of  California  [7], 

Based  on  these  observ  ations,  the  discharge  apparatus  has  been  specifically  configured  for  reac¬ 
tion  ion  etching  experiments  in  which  wafer  etching  stations  and  electrostatic  analyzers  are  located 
on  both  the  central  and  wall  electrodes.  The  immediate  objective  in  the  experiments  that  are 
presently  underway  is  to  examine  the  ion  energy  distributions  that  are  produced  in  CF4— 02  and 
('H4— H;  etching  gases,  their  effeclivness  in  etching  Si  and  SiO?  respectively,  and  their  influence  on 
substrate  surface  damage  as  measured  by  capacitance-voltage  methods. 
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SCIENTIFIC  OBJECTIVE: 

Excited  states  and/or  free  radicals,  produced  by  a  gas  discharge  or  an  intense  laser  beam,  play 
a  critical  role  in  the  etching  and  deposition  of  semiconductor  material.  One  case,  dry  processing 
with  plasmas,  has  become  the  "standard"  industrial  tool,  and  much  of  the  imaginations  involving 
laser  interactions  are  beginning  to  bear  fruit.  Inasmuch  as  these  sources  of  radicals  and  excited 
states  produce  an  environment  which  is  far  from  that  described  by  equilibrium  thermodynamics, 
one  can  expect  unique,  puzzling,  and  yet  most  interesting  processes  resulting  from  the  use  of  these 
techniques.  Our  goals  have  been  and  continue  to  be  the  understanding  of  the  complex  processes 
responsible  for  those  actions  and  a  development  of  a  scientific  model  which  quantifies  the  interac¬ 
tion  between  the  source  (electrons  or  photons),  the  donor  gas.  and  the  semiconductor  surface.  We 
will  continue  the  basic  study  of  the  kinetic  processes  occurring  in  these  environments,  but  we  also 
plan  to  combine  the  two  technologies. 

A  further  objective  is  to  exploit  novel  photochemical  processes  in  the  gas  phase  for  the  growth 
of  metal  and  semiconductor  films.  Particular  emphasis  will  be  placed  on  the  use  of  multiphoton 
ionization  (\1PI)  of  Column  IIIA  metal  alkyls  for  the  growth  of  compositionally  and  doping- 
modulated  III- V  film  structures. 

A  specific  target  for  this  next  period  is  to  combine  the  discharge,  with  its  ability  to  produce 
copious  quantities  of  free  radicals  in  almost  anything  rather  easily,  with  the  laser,  which  has  the 
advantage  of  energy  selectivity  and  directionality,  to  gain  further  insight  into  these  complex 
processes. 


SUMMARY  OF  RESEARCH: 

Plasma  Processing 

Completed  Work  --  During  the  past  year  we  have  completed  our  study  of  the  effects  of  an 
Oxygen  additive  on  the  silicon  etching  characteristics  of  discharges  in  \F3.  This  work  was  pub¬ 
lished  in  the  J.  Electrochemical  Society.  While  there  are  still  unknowns  about  the  chemistry  of 
such  discharges,  it  appeared  to  be  the  optimum  time  to  concentrate  our  efforts  on  the  effects  of  RF 
modulation  on  the  growth  of  a-Si:H. 


Current  work  —  The  past  year's  work  has  documented  the  parameter  space  for  the  growth  of 
hydrogenated  low  bandgap  amorphous  silicon  films  using  a  modulated  RF  discharge.  A  side-bv- 
side  comparison  of  film  properties  of  those  grown  with  conventional  C W  RF  discharges  and  those 
with  a  modulated  one  shows  the  superiority  of  the  latter  films:  lower  bandgap.  higher  density, 
higher  activation  energy,  and  as  good,  if  not  better,  photoconductivity. 

The  striking  difference  between  these  films  is  best  illustrated  by  Figure  1.  which  is  a  Tauc 
Plot  showing  the  optical  absorption  as  a  function  of  the  photon  energy.  The  extrapolation  of  the 
straight-line  portion  back  to  the  intercept  is  the  bandgap  and  thus  there  is  a  significant  reduction  in 
it  by  the  simple  expedient  of  square-wave-modulating  the  RF  discharge  excitation.  This  is  the  only 
change:  everything  else  is  the  same:  system-parallel  plane  excitation:  pressure  *-l  torr.  gas  mix¬ 
ture  -  ().5rr  SiH4  in  He:  flow,  peak-to-peak  RF  voltage;  peak  power:  and  substrate  temperature. 

Table  1  presents  this  side-by-side  comparison  of  the  films  grown  with  a  modulated  discharge 
to  that  using  a  CW  one  with  the  major  variable  between  the  entries  being  the  substrate  tempera¬ 
ture  One  usually  requires  a  high  substrate  temperature  ~250°C  in  order  to  reduce  the  bandgap— a 
trend  which  is  observable  from  the  table.  But  also  evident  is  the  fact  that  the  bandgap  is  insensi¬ 
tive  to  substrate  temperature  for  modulated  discharges.  Indeed,  the  films  grown  at  room  tempera¬ 
ture  in  this  manner  are  just  about  as  good  as  those  with  T  =  250°C.  This  could  be  technologically 
important  because  it  would  permit  a  nearly  infinite  variety  of  substrates  for  practical  applications. 

Unfortunately,  our  scientific  base  for  the  effect  is  still  in  its  formative  phase.  We  have  docu¬ 
mented  some  very  important  facts  about  such  discharges,  which  have  led  to  reasonable,  but.  as  yet. 
unproven  hypotheses.  For  instance,  modulation  appears  to  change  the  kinetics  somewhat  of  any 
plasma  processing  step,  etching  or  growth.  This  suggests  that  ion  bombardment  may  be  preparing 
the  surface  sites  during  the  "on"  time  and  the  surface  chemistry  dominating  during  the  "off"  period. 

We  do  know  that  a  copious  amount  of  disilane  is  produced  in  these  discharges  and  that  it  is 
not  responsible  for  the  very  contorted  time  history  of  the  bulk  electron  density.  We  have  also 
made  some  detailed  measurements  on  the  electron  kinetics  of  such  discharges.  Our  major  task  for 
the  coming  year  is  to  mold  these  separate  bits  of  information  into  a  reasonable  model. 

Laser  Growth  of  Films 

In  the  past  year,  the  experimental  effort  has  been  two-pronged.  Laser  photochemical  vapor 
deposition  (LPVD)  experiments  have  continued  in  which  Si2Hb  or  GeH4  is  pyrolvzed  while  the  sub¬ 
strate  is  illuminated  with  fluence  (  <  10  mJ  —  cm'3)  of  ultraviolet  laser  radiation.  The  film  growth 
rate  is  clearly  enhanced,  even  for  very  low  optical  fluences  (1—2  mJ  —  cm-2).  TEM  and  S\MS  stu¬ 
dies  are  being  carried  out  of  the  structural  (crystalline)  and  chemical  characteristics  of  the  grown 
films.  By  altering  the  laser  wavelength,  repetition  frequency,  and  pulse  energy,  it  has  also  been 
demonstrated  that  the  species  primarily  responsible  for  film  growth  can  be  determined  solely  by 
the  laser,  and  the  reactor  chemistry  changed  at  will. 

In  the  second  set  of  experiments,  a  turbomolecular  pumped  reactor  suitable  for  the  growth  of 
metal  or  semiconductor  films  from  atomic  ions  produced  by  multiphoton  ionization  has  been  con¬ 
structed.  We  are  currently  investigating  the  scaling  of  the  A1  ion  production  rate  with  laser 
intensity  when  ALlCHj)^  is  ionized  in  the  blue.  Similar  studies  will  soon  be  underway  with  gal¬ 
lium. 


y . <••••••: . € . . 


10  . 

8  . 

-  6  . 


. i . pqmiJUfe 

!  I  i  /  /  cw! 


2  3 

photon  energy  (eV) 


Fig  1  A  Tauc  plot  showing  the  typical  differences  between  films  deposited  from  continuous 
(CW)  and  modulated  (SQW'M)  discharges  (200  Hz  modulation  frequency) 


Table  1.  The  optical  and  electrical  properties  of  films  produced  in  CW  and  SQWM  discharges  (200 
Hz  modulation  frequency).  (*  data  unavailable) 
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SCIENTIFIC  OBJECTIVE: 

Our  principal  goal  on  the  investigation  of  monolithic  microwave  integrated  circuit  (MM1C)  is 
to  integrate  active  solid  state  devices  with  microstrip  antennas.  During  the  past  year,  progress  has 
been  made  in  the  study  of  propagation  characteristics  of  guidance  structures,  the  coupling  between 
waveguides  and  the  integration  of  solid-state  amplifiers  with  microstrip  antennas. 

SUMMARY  OF  RESEARCH: 

Propagation  Characteristics  of  Guidance  Structures 

The  theoretical  modeling  in  obtaining  the  propagation  characteristics  for  a  guidance  structure 
would  permit  the  design  of  a  proper  matching  network  to  the  device  and  optimum  feeding  to  the 
antenna.  Since  the  analysis  of  the  guidance  problem  plays  an  important  role  in  the  design,  consid¬ 
erable  amount  of  effort  must  be  expended  on  the  numerical  solution.  A  computer  program  using 
the  finite  element  method  has  been  developed  to  analyze  the  propagation  constants  for  dielectric 
waveguide  and  microstrip  lines.  The  program  has  the  advantage  of  being  able  to  handle  compli¬ 
cated  guidance  geometries  as  compared  to  conventional  techniques  such  as  the  effective  index 
method,  and  it  gives  better  results  near  cutoff  frequency.  However,  a  large  amount  of  computa¬ 
tional  time  and  memory  space  is  required  when  executing  the  program.  Ten  hours  of  supercom¬ 
puter  CPU  time  have  been  granted  to  work  on  this  project,  and  future  analyses  are  expected  to  be 
done  on  the  CRAY. 

Coupling  between  Waveguides 

During  the  past  year,  we  have  been  studying  the  applications  of  the  strongly  coupled-mode 
theory  to  optical  waveguide  couplers.  The  field  distribution  and  the  propagation  constants  were 
evaluated  using  the  coupled-mode  theory  and  compared  with  those  obtained  from  the  numerical 
approach  and  the  weakly  coupled-mode  theory.  The  differences  between  these  data  were  investi¬ 
gated.  and  we  conclude  that  the  strongly  coupled-mode  theory  is  a  better  method  as  compared  to 
the  weakly  coupled-mode  theory  under  either  strong  or  weak  coupling  conditions.  The  theory  was 
excellent  in  calculating  the  wave  distribution  and  the  energy  levels  for  coupled  quantum  wells 
under  equilibrium  conditions  and  under  the  influence  of  an  applied  electric  field.  The  coupled 
quantum  well  model  is  useful  for  novel  electronic  devices  such  as  resonant  tunneling  diodes,  which 
have  the  potential  for  millimeter-wave  generation. 
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Integration  of  Amplifier  to  Microstrip  Antennas 

Our  current  research  invokes  the  use  of  the  coplanar  transmission  line  as  a  means  of  monol- 
ilhically  integrating  solid-stale  dec  ices,  such  as  phase  shifters  and  I  L  IN.  ccith  an  array  ol  micros- 
trip  patch  antennas.  Due  to  the  novelty  of  the  coplanar  line  for  millimeter-cvave  applications, 
much  of  our  efforts  have  been  spent  on  analytically  and  experimentally  characterizing  coplanar 
structures  for  amplifier  design  as  well  as  lor  electromagnetically  feeding  the  antenna  array. 

We  have  dec  eloped  a  computer  program  cc  hich  can  quasistatically  calculate  the  impedance 
and  eileclice  dielectric  constant  ol  buried  and  unburied  coplanar  and  microstrip  transmission  lines. 
We  are  currently  extending  our  theoretical  analysis  to  a  full-wave  analysis  so  that  the  frequency 
dependence  ol  the  lines  may  be  determined,  (Experimentally  .  eve  have  demonstrated  the  feasibility 
ol  using  a  coplanar  line  in  monolithic  design  by  fabricating  and  testing  an  integrated 
amplifier  patch  antenna  prototype  using  coplanar  lines.  Currently,  eve  are  making  a  comparative 
studs  o!  several  different  methods  of  feeding  microstrip  antenna  arrays  using  microstrip  and 
coplanar  slot  line  feeds.  This  includes  comparisons  of  cross  polarization,  directivity,  input 
impedance  and  parasitic  radiation  from  the  coplanar  line  and  slot  line  In  addition,  we  are  develop¬ 
ing  a  novel  measurement  technique  which  can  be  used  to  characterize  accurate  coplanar  fixtures. 
Such  characterization  is  essential  to  the  accurate  measurement  of  transistor  S-parameters.  antenna 
input  impedance,  and  coplanar  discontinuities  lor  amplifier  design  This  measurement  technique 
also  has  potential  significance  in  the  field  of  microwave  measurements  due  to  its  simplicity  and 
cost 
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SCIENTIFIC  OBJECTIVE: 

There  are  three  principal  objectives  of  this  effort.  The  first  is  to  develop  efficient  integral 
equation  and  differential  equation  techniques  for  solving  the  problems  of  electromagnetic  scatter¬ 
ing.  and  coupling  into  complex  structures.  The  second  is  to  study  techniques  for  reducing  the  radar 
cross-section  of  targets  of  different  shapes.  The  third  is  to  analyze  frequency  selective  surfaces 
(FSS)  for  radomes  and  other  antenna  applications. 

SUMMARY  OF  RESEARCH: 

During  the  past  twelve  months,  we  have  successfully  solved  a  number  of  radar  scattering 
problems  for  bodies  whose  dimensions  are  large  compared  to  the  wavelength  of  the  illuminating 
field.  We  have  investigated  both  the  iterative  methods,  e  g.,  the  conjugate  gradient  method,  the 
spectral  iterative  technique  and  the  generalized  Neumann  series,  and  noniterative  techniques,  e.g.. 
the  spectral  Galerktn  method  and  the  moment  method.  We  have  carried  out  a  comparative  study  of 
the  scope  of  application  and  the  numerical  efficiency  of  these  techniques  by  solving  some  represen¬ 
tative  test  problems.  Two  examples  of  such  test  problems  are  a  large  thin  plate  and  a  long  missile¬ 
shaped  body  of  revolution.  An  important  aspect  of  the  RCS  study  has  involved  the  investigation  of 
the  multiple  incident  angle  problem  where  the  aspect  angle  of  the  incident  fieid  is  varied  over  a 
wide  angular  range,  e.g..  from  normal  incidence  to  grazing.  Because  the  body  dimensions  can  be 
large  compared  to  the  wavelength,  the  number  of  unknowns  employed  in  the  representation  of  the 
induced  current  on  the  scatterer  can  reach  up  to  several  thousands  when  subdomain  basis  functions 
are  employed.  The  resulting  equations  are  typically  solved  using  one  of  the  iterative  techniques 
because  their  use  can  circumvent  the  storage  problem  usually  associated  with  conventional  matrix 
elimination  schemes  when  a  large  number  of  unknowns  are  involved. 
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Despite  the  storage  advantage,  the  iterative  techniques  do  suffer  from  one  serious  drawback  in 
that  they  are  very  inefficient  for  solving  the  multiple  excitation  problem.  For  this  reason,  we  have  j 

also  investigated  the  possibility  of  reducing  the  number  of  unknowns  in  the  representation  of  the 
current  by  employing  the  entire  domain  type  of  basis  functions.  We  have  demonstrated  that  a 
reduction  of  the  number  of  unknowns  by  an  order  of  magnitude,  or  even  more,  is  possible  by  using 
an  appropriate  set  of  entire  domain  functions.  Using  these  functions,  we  have  been  able  to  solve 
plate  sizes  as  large  as  six  wavelength  square,  with  less  than  only  300  unknowns,  i.e..  with  the 
number  of  unknowns  at  least  an  order  of  magnitude  less  than  that  required  in  the  subdomain  for¬ 
mulation  We  are  currently  studying  the  susceptibility  of  the  solutions  generated  in  this  manner  to 
the  computer  round-off  errors.  This  study  is  important  because  the  condition  numbers  of  matrices 
generated  bv  using  the  entire  domain  basis  functions  is  often  much  larger  than  those  obtained  when 
the  subdomain  basis  functions  are  employed.  Hence,  the  matrix  for  the  entire  domain  case  is  usu¬ 
ally  ill-conditioned.  Methods  for  alleviating  this  difficulty  are  currently  under  investigation  and 
the  preliminary  results  look  very  promising. 

In  addition  to  the  plate  problem,  we  have  been  investigating  the  body  of  revolution  (BOR) 
problem,  which  plays  an  important  role  in  modeling  a  number  of  canonical  radar  targets.  We  have 
developed  the  theoretical  analysis  and  the  computer  program  based  on  this  analysis  for  a  BOR 
coated  with  one  or  more  layers  of  lossy  dielectrics.  Numerical  results  have  been  generated  for 
differentlv  shaped  bodies,  e.g..  cylinders,  spheres,  cones,  and  combinations  thereof,  and  the  results 
have  been  compared  with  those  reported  elsewhere.  Body  sizes  up  to  10  wavelengths  have  been 
analyzed  with  good  success. 

The  problem  of  mutual  coupling  between  antennas  mounted  on  arbitrarily  shaped  cylinders  is 
an  important  one  from  the  point  of  view  of  estimating  the  degree  of  interference  between  two 
antennas  located,  say.  on  the  fuselage  or  the  wing  of  an  airplane.  This  problem  can  be  difficult  to 
solve  using  the  asymptotic  techniques,  e.g.,  the  geometrical  theory  of  diffraction  (GTD).  even  when 
the  size  of  the  cylinder  is  large  compared  to  the  wavelength.  We  have  developed  numerical  tech¬ 
niques  for  solving  this  problem  and  have  generated  results  for  some  representative  choice  of 
parameters  for  the  antennas  and  the  cylinder. 

Finally,  we  have  developed  iterative  techniques  for  solving  the  problem  of  frequency  selective 
screens  comprised  of  a  periodic  array  of  arbitrarily  shaped  patches  or  apertures.  The  program  is 
capable  of  handling  a  large  number  of  unknowns  and  FSS  elements  with  arbitrary  geometries. 

W'e  are  also  studying  the  effect  of  truncating  the  doubly  periodic  screens  and  are  evaluating 
the  effects  of  this  truncation  on  the  scattering  properties  of  this  screen.  Another  aspect  of  our  study 
involves  the  investigation  of  the  out-of-band  scattering  characteristics  of  these  screens.  These  j 

characteristics  are  of  interest  from  the  point  of  view  of  designing  satellite  antenna  systems  with  1 

low  RCSs.  We  are  also  investigating  the  effect  of  curvature  on  the  frequency  selective  properties 
and  the  scattering  characteristics  of  FSS  screens.  The  curved  FSS  problem  is  important  in  the  design  , 

of  radomes.  .  j 

A  comparative  study  of  various  techniques  for  solving  large-bodv  scattering  problems  has 
been  made,  and  some  definitive  conclusions  regarding  the  relative  advantages  and  disadvantages  of 
these  techniques  have  been  reached.  A  spectral  Galerkin  approach  has  been  developed  for  the  pur-  i 

pose  of  reducing  the  number  of  unknowns  and  has  been  successfully  applied  to  large-body¬ 
scattering  problems.  The  problem  of  estimating  the  mutual  coupling  between  antennas  mounted  on 
airplanes  or  missiles  has  been  studied.  Frequency  selective  surfaces  with  arbitrary  elements  have  ‘ 

been  investigated  and  the  study  of  truncated  as  well  as  conformal  FSSs  has  been  initiated. 
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SCIENTIFIC  OBJECTIVE: 

This  unit  seeks  to  develop,  model,  and  analyze  efficient,  high-performance,  and  testable  com¬ 
puter  architectures  which  will  exploit  compound  semiconductor  technology.  There  has  been 
increasing  interest  in  the  use  ol  compound  semiconductor  technologies,  such  as  gallium  arsenide 
(GaAs).  as  a  vehicle  to  develop  high-performance  electronic  systems.  Various  studies  have  indi¬ 
cated  that  GaAs  devices  will  become  a  larger  share  of  the  semiconductor  market  in  the  coming 
years.  In  addition  to  the  high  performance  possible,  the  capability  for  greater  radiation  hardness  in 
the  compound  semiconductor  technology  is  of  great  interest  to  military  systems.  We  have 
identified  several  aspects  of  computer  architecture  for  particular  attention,  due  both  to  their  emerg¬ 
ing  importance  from  a  technology-driven  point  of  view,  and  to  the  lack  of  known  structures  or 
analysis  techniques  for  meeting  our  objective. 


SUMMARY  OF  RESEARCH: 

Memory  Organizations 

Trace-driven  simulation  is  a  simple  way  of  evaluating  cache  memory  systems  with  varying 
hardware  parameters.  To  evaluate  realistic  workloads,  simulating  even  a  few  millions  of  addresses 
is  not  adequate  and  such  large-scale  simulation  is  impractical  from  the  consideration  of  space  and 
time  requirements.  Therefore,  simply  reporting  some  more  simulation  results  on  cache  memories, 
which  is  found  in  profusion  in  existing  literature,  is  not  the  goal  of  this  research.  The  goal  of  this 
research  is  to  develop  new  methods  of  simulation  based  on  statistical  techniques  for  decreasing  the 
need  for  large  trace  measurements  and  to  predict  true  program  behavior.  In  the  proposed  method, 
sampling  techniques  are  applied  while  collecting  the  address  trace  from  a  workload.  The  cost  of  the 
sampling-based  simulation  method  in  terms  of  both  memory  space  and  computer  time  is 
significantly  lower  than  conventional  techniques.  For  a  trace  length  of  5  million  references,  this 
method  typically  requires  collecting  only  7%  or  less  of  the  trace  for  smaller  caches.  For  large 
caches,  much  longer  segments  of  traces  are  required  by  conventional  methods.  The  new  method 
estimates  the  true  value  of  miss  ratio  with  much  lower  overhead.  This  technique  gives  accurate 
estimates  of  both  the  mean  value  and  the  distribution  of  the  miss  ratio. 

•  The  basic  sampling  technique  that  has  been  established  is  applicable  to  small  caches.  It  is 
based  on  the  assumption  that  caches  are  Hushed  totally  during  context  switches.  This  is  essentially 
true  for  smaller  caches.  The  distribution  of  the  miss  ratio  of  the  samples  gives  a  reasonably 


accurate  estimate  of  that  of  the  cold-start  miss  ratio  over  the  entire  trace  for  the  same  context 
switch  interval  as  the  sample  size.  It  has  been  experimentally  shown  for  a  varied  class  of  programs 
that  about  35  samples  are  adequate  to  well  characterize  the  program  ixehav  lor  The  distribution  ol 
the  miss  ratio  obtained  by  the  sampling  technique  is  compared  with  that  trom  the  simulation  oi  the 
corresponding  continuous  trace  by  means  ol  statistical  techniques  to  validate  the  sampling-based 
simulation  method. 

The  method,  mentioned  above,  predicts  the  miss  ratio  oi  programs  with  context  switch  inter¬ 
val,  equal  to  the  sample  size.  A  model  is  proposer!  to  statistically  project  the  results  to  different 
context-sw  itch  intervals  1  rom  only  one  simulation  of  samples  ol  a  fixed  size.  The  validity  ol  this 
projection  has  been  experimentally  established  This  model  is  also  extender!  to  predict  t lie  miss 
ratio  of  continuous  tasks  when  the  program  runs  to  completion  without  any  context  switching. 

The  large  caches  are  becoming  very  common  and  their  behavior  can  only  be  studied  by  using 
massively  long  traces.  But.  the  above  sampling-based  methods  are  not  applicable  to  large  caches 
where  significant  amounts  of  data  are  retained  across  the  context  switches  Therefore,  a  new  con¬ 
cept  of  primed  cache  is  introduced  to  simulate  large  caches  by  the  sampling-based  method  1  ms  is 
achieved  by  simulating  only  those  portions  of  the  cache  whose  true  state  can  be  recreated  even 
from  the  sampled  trace  Both  direct-mapped  and  set-associative  caches,  as  large  as  128  Kbytes, 
have  been  studied  using  these  methods. 

Parallel  Pipelined  Architectures 

We  have  been  studying  new  techniques  for  accurate  performance  estimation  of  VLSI  chips. 
This  is  a  crucial  problem  in  the  design  of  very  high-speed  systems  which  include  parallel  pipelined 
architectures  Current  techniques  are  limited  in  two  respects.  First.  SPICF-tvpe  simulations  are 
impractical  due  to  the  excessive  CPU  time  required  to  simulate  large  chips.  Second  less  accurate 
timing  simulators  tend  to  calculate  the  absolute  worst-case  path,  wothoui  regard  to  whether  this 
path  is  ever  actually  energized  We  are  developing  a  hierarchical  timing  simulation  system  that  w  ill 
use  high-level  information  (at  the  register-transfer  level,  for  example)  to  determine  the  worst-case 
paths  for  the  defined  operations  of  the  system.  The  values  of  the  worst-case  delays  will  be  deter¬ 
mined  using  high-level  floorplan  information  and  intra-module  timing  information.  The  layout 
information  will  allow  the  simulation  system  to  use  transmission  line  models  to  estimate  commun¬ 
ication  delays  between  modules  The  use  of  transmission  line  models  is  necessary  since  the 
increased  circuit  switcning  speeds  of  GaAs  and  other  technologies  result  in  long  interconnections 
behaving  like  transmission  lines. 
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SCIENTIFIC  OBJECTIVE: 

This  unit  is  concerned  with  exploring  new  directions  in  fault  tolerant  computing  that  will  be 
useful  in  understanding  the  basic  principles  in  design,  testing,  error  prediction  and  circumvention 
for  reliable  complex  computer  systems.  Basic  research  is  being  performed  in  the  three  areas  of 
parallel  test  pattern  generation,  reliable  distributed  database  systems  for  real-time  applications, 
and  intelligent  error  prediction  and  circumvention. 


SUMMARY  OF  RESEARCH: 

Parallel  Algorithms  for  Automatic  Test  Generation 

Our  initial  investigations  into  the  design  of  parallel  algorithms  for  automatic  test  pattern  gen¬ 
eration  on  multiprocessor  systems  have  shown  a  lot  of  promise.  We  have  identified  some  tech¬ 
niques  for  possible  parallelization  which  we  are  going  to  implement  and  evaluate  in  the  future.  ( 1 ) 
Partitioning  a  large  circuit  into  blocks  which  are  assigned  to  different  processors.  (2)  Performing 
forward  implications  in  parallel  by  guessing  signal  values  on  lines  based  on  heuristics;  the  correct¬ 
ness  of  the  guessed  values  are  checked  by  interprocessor  messages.  (3)  Using  multiple  backtrace 
procedures  to  trace  signal  determination  backwards  on  multiple  paths.  (4)  Performing  D-drive 
through  all  the  gates  of  D-frontier  simultaneously. 

Reliable  Distributed  Database  Systems 

Research  concerning  concurrent  detection  and  recovery  of  structure  and  semantic  errors  in 
database  systems  has  resulted  in  a  virtual  backpointer  approach  for  the  on-line  detection  of  errors 
in  linked  data  structures.  Both  Virtual  Double-Linked  Lists  and  B-trees  with  Virtual  Backpointers 
have  been  developed  based  on  the  virtual  backpointer  concept.  The  data  structures  are  appropriate 
for  real-time  applications  since  they  allow  concurrent  error  detection  in  0(1)  time  with  little  or  no 
storage  overhead  and  few  extra  node  accesses.  For  purposes  of  analysis,  checking  window  and  lock 
and  key  concepts  have  also  been  introduced.  These  techniques  have  been  implemented  in  a  data 
structure  software  package  and  are  currently  undergoing  analysis.  The  implementations  are  being 
investigated  for  performance  as  well  as  reliability  attributes  in  both  shared  memory  multiprocessor 
and  distributed  work-station  environments. 


Intelligent  Error  Prediction  and  Circumvention 

A  methodology  of  model  development  from  the  raw  error-data  to  the  estimation  of  cumula¬ 
tive  reward  has  been  developed  The  results  provide  an  understanding  of  the  different  types  ot 
errors  and  recovery  processes  The  measured  data  show  that  the  holding  times  in  key  operational 
and  failure  states  are  not  simple  exponential  anil  that  a  semi -Markov  process  is  necessary  to  model 
the  system  behavior  A  sensitivity  analysis  has  also  been  performed  to  investigate  the  significance 
of  using  a  semi-Markov  process  (as  opposed  to  a  Marko^  process)  to  model  the  measured  system. 
Results  show  that  the  Markov  process  frequently  overestimates  the  error  probabilities  and 
underestimates  the  second  moment  ot  the  first  passage  lime  to  the  error  states. 

Based  on  the  rnodei  methodology  .  a  measurement -based  perl ormabtlily  model  has  been 
developed.  A  reward  rate  is  used,  which  is  defined  as  a  function  ot  service  rate  and  the  error  rate 
in  each  stale,  to  estimate  the  perlormability  ot  the  sy  stem  and.  to  depict  the  cost  of  different  error 
types  and  recovery  procedures.  Results  show  that  not  recovering  from  software  errors  degrades 
system  effectiveness  considerably  compared  with  the  effect  of  not  recovering  from  disk  errors  if  we 
assume  the  system  continues  to  provide  service  in  a  degraded  mode.  However,  if  we  assume  that 
both  software  and  disk  errors,  result  in  system  failed  then  the  effect  due  to  the  disk  errors  is  more 
severe  than  that  due  to  the  software  errors.  A  software  reliability  model  has  also  been  built  based 
on  the  low-level  error  data  from  the  MVS  operating  system  running  on  an  IBM  3081  machine.  The 
semi-Markov  model  developed  provides  a  quantification  of  system  failure  characteristics  and  the 
interaction  between  different  types  of  errors. 
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SCIENTIFIC  OBJECTIVE: 

The  analysis  and  the  design  of  efficient  computation  techniques  has  been  for  several  years  one 
of  the  most  vigorous  areas  of  research  in  information  processing  Its  impact  has  been  not  only  on 
the  availability  of  better  computational  methods  for  a  number  of  significant  applications,  but  also 
on  the  awareness  in  the  computing  community  of  the  crucial  importance  of  algorithmic  design.  By 
investigating  both  upper  and  lower  bounds  to  the  resource  requirements  of  specific  applications, 
this  discipline  effectively  develops  a  methodology  aiming  at  a  quantitative  formulation  of  the 
inherent  difficulty  of  problems  in  the  appropriate  computation  models. 

The  involvement  of  cur  research  group  with  this  topic  is  well  established  and  dates  back  to 
the  early  seventies.  Our  record  over  this  period  shows  that  our  focus  (as  well  as  that  of  our  peer 
community)  has  been  adjusting  dynamically,  as  changes  in  technologies  have  modified  the  general 
research  horizon  and  the  perception  of  relevance.  The  most  significant  change  of  this  type  has  been 
the  advent  of  Very-Large-Scale-Integration  (VLSI),  which  has  profoundly  influenced  essentially 
every  facet  of  our  current  research  interests.  The  advent  of  VLSI  is  important  in  two  major 
respects:  one  is  the  present  possibility  to  realize  massively  parallel  computers,  the  other  is  the 
introduction  of  criteria  of  complexity  (the  VLSI  "model")  which  take  into  account  the  design  rules 
dictated  by  the  new  technology.  Therefore,  while  maintaining  our  interest  for  problems  in  the 
more  traditional  areas  (Von  Neumann  computations),  our  current  work  emphasizes  research  on  the 
potentially  revolutionary  domain  of  concurrent  computation 


SUMMARY  OF  RESEARCH: 

Parallel  Computation  and  Graph  Algorithms 

In  parallel  computation,  we  have  been  looking  at  fast  parallel  algorithms  for  combinatorial 
problems.  We  have  developed  a  fast  parallel  algorithm  for  the  linear  programming  problem  with 
two  variables  per  inequality  and  two  variables  in  the  optimization  function.  The  parallel  algorithm 
runs  in  polylog  time  with  a  subexponential  number  of  processors. 
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In  continuation  of  our  work  on  graph  algorithms:,  we  have  designed  several  parallel  XC  algo¬ 
rithms  for  reducible  flow  graphs.  We  have  XC  algorithms  lor  recognizing  rfg's.  lor  finding  a  depth 
first  search  tree  in  an  rfa  anti  for  finding  a  minimum  feedback  vertex  set.  We  have  also  developed 
an  RXC  algorithm  for  finding  a  minimum  weight  feedback  arc  set  in  arc-weighted  rfg’s. 

Vertex  connectiv  uy  is  an  important  graph  property.  We  have  developed  a  new  algorithm  for 
testing  vertex  four-connectivity  and  for  finding  all  separation  triplets  in  a  triconnected  graph  in 
( )(  n  ’ 1  sequential  time,  anti  in  ( )(  log’n )  parallel  time  on  a  CRCW  PRAM  with  O(n’)  processors.  This 
improves  prev  ious  bounds  lor  the  problem  for  both  the  sequen'  al  and  parallel  cases. 

We  have  looked  into  efficient  parallel  circuits  and  algorithms  for  the  problem  of  finding  the 
n-bit  result  of  d:  •  iding  one  n-bit  number  by  another.  We  have  developed  a  logspace  unilorm  lam- 
:1c  ol  circuits  for  dr.  is  ion  of  depth  D(  n  )=()(,(  l/5-*)lognloglogn)  and  size  S(  n  )=( )( (  1 /S'j  n1 "  ).  lor 
any  5  >  O.  This  translates  into  a  uniform  parallel  algorithm  on  a  shared  memory  machine  (PRAM) 
with  bit  operations  and  exclusive  memory  writes  with  parallel  time  D(n)  using  0(S(n))  processors. 
It  also  translates  into  a  parallel  algorithm  for  a  concurrent  write  PRAM  with  parallel  time 
C)(D(  n)/loglogn  )  using  ()( St  n ) )  processors.  We  have  also  obtained  a  polynomial-time  unilorm  fam¬ 
ily  of  circuits  for  div  ision  of  depth  ()(( l/S2)  logn)  and  size  0(S(n)).  Finally  we  have  looked  to 
small  depth  circuits  for  the  evaluation  of  arithmetic  expressions.  We  have  some  preliminary  results 
on  this  problem. 

Models  of  Parallel  Computation 

The  most  important  model  of  parallel  computation  is  the  Parallel  Random  Access  Machine 
(PRAM),  which  comprises  multiple  processors  that  execute  instructions  synchronously  and  share  a 
common  memory.  This  model  provides  the  foundation  for  the  design  of  highly  parallel  algorithms. 
Because  of  the  widespread  use  of  this  model,  an  understanding  of  its  capabilities  is  crucial.  The 
usual  PRAM  model  has  the  following  instructions:  bitwise  boolean  operations:  integer  addition, 
subtraction,  and  comparisons  with  conditional  jumps:  indirect  read  and  write.  We  have  studied  the 
computational  capabilities  of  deterministic  PRAMs  with  additional  instructions. 

We  proved  that  polynomial  space  on  Turing  machines  ( PSPACE )  is  equivalent  to  polynomial 
time  on  PRAM s  with  unit-time  multiplication  and  on  PRAM s  with  unit-time  unrestricted  left 
shift  and  right  shift.  These  results  are  surprising  for  two  reasons.  First,  for  the  sequential  random 
access  machine  running  in  polynomial  time,  adding  unit-time  multiplication  or  unit-time  unres¬ 
tricted  left  shift  increases  its  power  from  PTIME  to  PSPACE.  whereas  adding  one  of  these  opera¬ 
tions  to  a  PRAM  does  not  increase  its  power  beyond  a  polynomial  amount.  Second,  despite  the 
potential  speed  offered  by  massive  parallelism,  a  PRAM  with  unit-cost  multiplication  or  unres¬ 
tricted  shifts  is  at  most  polynomiallv  faster  than  a  sequential  random  access  machine  with  the 
same  additional  operation. 

VLSI  Layout 

We  have  developed  new  algorithms  for  multilayer  channel  routing  and.  moreover,  we  now 
have  a  unified  framework  in  which  many  previously  known  results  can  be  obtained.  We  have  con¬ 
sidered  various  models,  including:  Manhattan,  knock-knee,  unit  vertical  overlap,  routing  with 
diagonal  wires,  and  unrestricted  overlap.  In  addition,  we  are  investigating  the  possibility  of  having 
"stacked  pins."  Some  lower-bound  results  have  also  been  obtained. 

We  have  also  looked  al  the  question  of  delay  in  VLSI  layouts  and  the  area  penalty  caused  by 
making  iocal  expansions  in  layouts.  We  have  obtained  tight  bounds  for  these  for  various  types  of 
layouts.  In  particular,  we  show  that  reducing  capacitive  delay  in  layouts  by  introducing  drivers 
can  sometimes  cause  a  large  area  increase. 

Complexity  of  Graph  Problems 

Recent  advances  in  graph  theorv  and  graph  algorithms  may  have  profound  implications  for 
the  field  ot  complexity  theory  Work  by  Robertson  anti  Seymour  has  provided  powerful  new  tools 


to  show  problems  in  P  by  proving  only  the  existence  of  polynomial  time  decision  algorithms.  We 
have  applied  these  techniques  to  show  a  variety  of  problems  to  be  polynomial  time  solvable;  some 
were  not  previously  known  even  to  be  decidable,  others  were  known  to  be  in  P.  but  low -order 
pol>  nomial  time  solutions  were  not  known.  Examples  come  from  such  diverse  areas  as  topological 
embeddings,  reliability  and  fault-tolerance,  and  V  LSI  applications,  as  well  as  straightforward 
graph  problems  Moreover,  a  number  of  the  problems  are  self-reducible:  answering  the  decision 
question  can  be  used  to  construct  a  solution. 
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SCIENTIFIC  OBJECTIVE: 

The  term  "sensor  array  system"  refers  to  a  large  class  of  remote  sensing  systems  in  which  data 
are  collected  and  recorded  by  many  independent  sensors,  or  by  one  sensor  that  is  moved  to  different 
spatial  positions.  The  recorded  data  are  processed  by  a  digital  array  algorithm  to  produce  a  high- 
resolution  object  function  Some  of  the  more  important  multi-sensor  array  systems  now  in  use  are 
synthetic  aperture  radar  (SAR).  computer-aided  tomography  (CAT),  and  beam-forming  sonar.  The 
objective  of  the  research  in  this  unit  is  to  develop  both  the  theory  and  computer  verification  of  new 
signal  processing  methods  with  the  goal  of  overcoming  current  limits  on  resolution  and  data 
throughput  rates  for  these  systems.  Specifically,  we  propose  to  develop  reconstruction  algorithms 
that:  (1)  can  achieve  high  resolution  from  a  limited  amount  of  sampled  data;  (2)  have  small 
enough  computational  complexities  that  real-time  implementation  is  feasible;  and  (3)  can  be  parti¬ 
tioned  into  sets  of  elementary  operations  suitable  for  VLSI  realization. 


SUMMARY  OF  RESEARCH: 

During  the  past  year  our  JSEP  research  has  concentrated  on  two  aspects  of  high-resolution 
sensor  array  processing:  adaptive  filtering  and  beamforming  algorithms,  and  reconstruction  algo¬ 
rithms  for  synthetic  aperture  radar  (SAR).  In  the  adaptive  area  we  have  continued  studying 
discrete  orthogonal  transforms  for  improving  convergence  properties  in  adaptive  finite  impulse 
response  (FIR)  digital  filters  which  operate  in  colored  noise  environments.  Also,  we  have 
transferred  some  of  our  results  to  the  problem  of  adaptive  beamforming  and  have  shown  that  the 
concept  oi  prewhitening  by  means  of  applying  an  orthogonal  transformation  can  improve  conver¬ 
gence  properties  of  adaptive  beamformers.  Related  to  SAR.  we  have  developed  a  fast  Hankel 
transform  algorithm  for  accurately  inverting  polar-format  SAR  data,  and  we  have  analyzed  the 
spatial  domain  effects  of  the  Fourier  domain  polar-to-Cartesian  interpolation  utilized  in  spotlight 
mode  SAR.  We  have  also  examined  the  possibilities  of  artificially  extending  (extrapolating)  the 
known  Fourier  data  for  improved  resolution  and  of  reconstructing  from  just  the  phase  of  the 
Fourier  data  to  simplify  computation. 


Adaptive  Filtering  and  Beamforming 

After  having  achieved  significant  improvements  in  FIR  digital  filters  with  the  Walsh- 
Hadamard  transform  (WHT)  [2],  we  designed  a  new  orthogonal  transform  called  the  Power-of- 
1'vvo  (P()2)  transform  that  has  a  similar  structure  to  the  WHT  hut  which  has  real-valued 
coefficients  that  are  simple  powers-of-2  [13].  Since  the  1*02  transform  is  a  generalization  of  the 
WHT.  it  was  expected  that  it  would  achieve  faster  convergence  rates  while  requiring  only  shifts 
and  adds  for  implementation:  this  expected  improvement  was  demonstrated  in  numerous  computer 
experiments  It  was  recently  shown  that  the  discrete  Hartley  transform  (DHT).  which  is  a  close 
relative  of  the  FFT.  preforms  belter  than  any  ol  the  FFT.  WHT.  or  P()2  transforms  and  is  perhaps 
the  strongest  candidate  for  real-time  adaptive  filtering  applications. 

In  broadband  beamforming,  each  sensor  in  the  array  is  typically  followed  by  an  FIR  digital 
filter  to  provide  the  necessary  frequency -dependent  array  weighting,  and  the  weights  can  be 
adjusted  adaptively  to  steer  the  nulls  of  the  beampattern  toward  any  desired  sources.  For  adaptive 
beamformers  based  on  the  l.MS  algorithm  in  a  multi-signal  environment,  convergence  is  data 
dependent  and  is  characterized  by  highly  disparate  modes,  resulting  in  slow  and  noisy  adaptation. 
Recently,  we  began  to  examine  the  feasibility  of  using  orthogonal  transforms  (discussed  above  for 
adaptive  filtering)  for  the  multichannel  wideband  beamforming  problem  to  improve  convergence 
properties  of  an  adaptive  beamformer.  A  special  class  of  beamformers.  the  generalized  sidelobe 
canceller,  was  studied  This  structure  converts  the  problem  of  linearly  constrained  adaptive  beam- 
forming  to  one  of  unconstrained  multiple  reference  noise  cancelling.  Applying  orthogonal 
transforms  to  this  multiple  input  case  is  nontrivial,  since  many  different  orderings  of  the  input 
data  vectors  are  possible  prior  to  transformation.  Overall,  it  was  demonstrated  by  computer  exper¬ 
iments  that  improved  convergence  rates  can  be  achieved  using  fixed  (time-invariant)  orthogonal 
transforms,  but  performance  is  highly  dependent  upon  input  data  and  the  array  processing  struc¬ 
ture  The  desired  result  of  identifying  a  single  transform  that  works  well  for  a  broad  class  of 
input  signals  and  different  numbers  of  sensors  has  not  yet  been  accomplished.  Further  research  is 
needed  to  develop  criteria  for  choosing  a  good  transform  from  among  the  many  possible  choices,  for 
ordering  the  data  vectors  that  are  presented  as  inputs  to  the  transform  algorithm,  and  for  evaluat¬ 
ing  the  results  on  an  analytical  basis. 

SAR  Reconstruction  Algorithms 

We  have  obtained  several  interesting  and  useful  results  on  algorithms  for  SAR  image  recon¬ 
struction.  In  spotlight  mode  SAR,  preliminary  processing  of  the  raw  data  provides  samples  of  the 
Fourier  transform  of  the  terrain  or  scene  reflectivity.  Unfortunately,  these  samples  are  available 
on  a  polar  grid  in  only  a  very  small  region  of  Fourier  space.  One  scheme  to  invert  the  Fourier  data 
is  to  first  interpolate  to  a  Cartesian  grid  and  to  then  apply  a  2-D  FFT.  In  [8]  we  studied  the  effect 
of  imperfect  Fourier  domain  interpolation  and  showed  that  serious  artifacts  may  arise  in  the  spatial 
domain.  Both  the  locations  and  amplitudes  of  these  artifacts  can  be  predicted.  In  [l]  we  considered 
the  possibility  of  extrapolating  the  Fourier  data,  to  increase  resolution,  prior  to  inversion.  It  was 
shown  that  conventional  extrapolation  algorithms  will  yield  increased  resolution  only  if  the  width 
of  the  known  Fourier  data  is  less  than  the  bandwidth  of  just  the  magnitude  of  the  scene 
reflectivity.  Publication  [3]  shows  that  for  scenes  of  high  contrast,  a  high  quality  SAR  image  can  be 
produced  from  just  the  phase  of  the  Fourier  data,  offering  the  possibility  of  considerably  reduced 
computation.  Finally,  in  [9]  a  new  algorithm  was  developed  for  efficiently  computing  a  series  of 
Hankei  transforms.  This  algorithm  was  subsequently  applied  in  a  Hankel  transform-based  algo¬ 
rithm  for  inverting  polar  Fourier  data  in  a  highly  accurate  manner.  Indications  are  that  the  Hankel 
transform  approach  is  very  competitive  in  terms  of  quality  and  computational  cost  with  the  more 
well-known  convolution-back-projection  algorithm. 
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SCIENTIFIC  OBJECTIVE: 

To  develop  high  throughput,  parallel  VLSI  structures  for  the  real-time  implementation  of 
high-resolution  algorithms  required  in  sensor  array  systems  such  as  synthetic  aperture  radars, 
beamforming  sonars,  direction  finders,  x-ray  cat  scanners,  and  acoustic  and  A'MR  imaging  systems. 
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SUMMARY  OF  RESEARCH: 

Research  on  this  new  JSEP  unit,  which  began  on  October  1.  1986.  has  concentrated  on  high- 
throughput  multiprocessor  VLSI  structures  for  some  common  linear,  shift-invariant  signal  process¬ 
ing  kernels  such  as  FIR  and  HR  digital  filtering,  on  systematic  design  methodologies  for  mapping 
signal  processing  algorithms  into  regular  and  modular  VLSI  structures,  and  on  studying  the  use  of 
residue  number  system  (RNS)  arithmetic  for  designing  a  fault-tolerant  arithmetic  cell  to  be  used  as 
a  standard  cell  in  the  design  of  VLSI  sensor  array  processors. 


i 
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We  have  designed  ultra-fast  structures  for  low-order  digital  filters,  using  a  large  number  of 
processors  to  generate  multiple  outputs  per  multiplication  cycle  (which  can  be  as  small  as  a  clock 
cycle).  These  structures  can  have  throughput  rates  many  times  higher  than  the  clock  frequency  on 
the  chip.  Earlier,  we  had  reported  work  on  modified  conventional  direct-form  FIR  and  HR  filter 
structures  for  high-throughput  VLSI  implementation  [ij.  Since  this  unit  began  in  October  1986.  we 
have  also  obtained  alternate  high-throughput  filter  structures  by  modifying  other  conventional 
filter  structures  such  as  the  cascade  and  parallel  forms.  It  is  expected  that  these  alternate  struc¬ 
tures  will  be  less  sensitive  to  finite-precision  errors.  These  new  structures  use  the  same  number  of 
processors  as  the  modified  direct  form  structures  of  [  1  ].  for  a  given  filter  order  and  a  given 
throughput  rate.  Conventional  parallel  and  cascade  forms  are  well  known  to  have  superior  finite- 
precision  behavior,  compared  to  the  conventional  direct  forms  ("conventional"  refers  to  structures 
that  produce  one  output  per  multiplication  cycle).  We  will  demonstrate  that  the  same  properties 
hold  for  the  new  high-throughput  structures.  Although  none  of  the  new  high-throughput  struc¬ 
tures  achieve  100%  processor  utilization  in  the  IIR  case,  they  all  achieve  100%  utilization  for  the 
FIR  application.  In  both  cases,  the  processor  utilization  factor  is  higher  than  that  achieved  by  the 
general  state-space  structures  that  have  received  wide  attention  in  the  literature. 


During  the  past  year  we  continued  our  study  of  systematic  design  methodologies  to  map  algo¬ 
rithms  into  regular  modular  structures.  We  have  studied  the  complexities  in  volume  and  wire 
length  of  mapping  problems  in  two  and  three  dimensions  [4.6],  studied  the  design  ot  macropipeline 
of  systolic  arrays  [<S],  and  investigated  the  mapping  of  dynamic  programming  algorithms  into  sys¬ 
tolic  arrays  [5],  We  have  also  continued  to  study  the  application  of  R\S  arithmetic  to  designing 
VLSI  modules  that  can  be  used  for  high-resolution  sensor  array  processing  [2,7.11.12].  A  semi- 
custom  integrated  circuit  module  was  designed  for  an  RMS  digital  tiller  of  order  X  using  the  IBM 
MVISA  ('AD  svstem  that  is  located  on  the  campus  ot  the  University  ot  Illinois  [  1 3].  (Mote  that 
this  work  was  done  largely  before  unit  17  began,  and  is  now  continuing  under  this  sponsorship.)  A 
programmable  4-bit  module  was  designed  using  finite  field  logarithm  addition  to  completely  elim¬ 
inate  multipliers  on  the  chip.  The  chip  design  was  done  with  MVISA.  which  then  was  used  to 
simulate  the  hardware  operation  and  estimate  physical  parameters.  MVISA  estimates  that  the  4  2 
mm  chip  consumes  89  mW  and  operates  at  a  system  cycle  frequency  of  10  mHz.  which  corresponds 
to  a  data-cycle  frequency  of  1.2  mHz.  The  design  used  529  out  of  a  possible  560  available  cells. 
The  chip  was  recently  fabricated  at  the  IBM  Manassas  facility  and  is  now  undergoing  testing  and 
verification. 
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SCIENTIFIC  OBJECTIVE: 

Our  group  has  been  and  is  in  the  forefront  of  research  on  a  wide  array  of  problems  in  adap¬ 
tive  systems:  some  of  the  problem  formulations,  methods  and  results  which  have  originated  in  our 
research  group  have  been  accepted  by  leading  researchers  worldwide  and  are  being  experimentally 
tested.  This  accumulated  experience  and  knowledge  in  the  field  puts  us  in  the  best  position  to 
address  the  underlying  fundamental  issues  and  to  develop  a  theory  for  adaptive  algorithms  which 
operate  in  partially  modeled  environments. 

Our  approach  embodies  the  philosophy  that  adaptive  systems,  be  they  adaptive  filters,  estima¬ 
tors  and  controllers,  or  diagnostic  and  tuning  software  for  automatic  quality  control  and  mainte¬ 
nance.  are  logical  on-line  extensions  of  well-tested  off-line  algorithms,  and  as  such,  should  inherit 
their  robustness  properties.  To  achieve  this  level  of  realism,  we  allow  our  convergence  and  stabil¬ 
ity  conditions  to  be  signal-dependent,  and  guarantee  them  in  the  domain  of  interest,  rather  than 
globally.  This  closes  the  gap  separating  the  realistic  conditions  from  the  mathematically  appealing 
global  results  derived  under  unrealistic  and  unverifiable  modeling  assumptions.  The  robustness 
properties  guaranteed  by  our  approach  will  compensate  for  the  loss  of  fictitious  globalitv. 

Toward  this  general  objective,  one  of  our  goals  in  this  research  project  is  to  investigate  the  key 
question  of  whether  there  are  self-tuning  algorithms  with  the  important  feature  that  the  algorithm 
will  automatically  tune  itself  to  a  correct  regulator  for  a  given  unknown  system.  Earlier  versions 
of  such  algorithms  have  had  a  big  impact  on  the  practical  control  of  a  certain  class  of  systems,  and 
are  also  beginning  to  be  fabricated  on  chips  and  installed  in  industrial  plants,  vehicles,  ships,  and 
other  communication  and  control  equipment. 

Our  second  goal  is  to  examine  the  behavior  of  the  algorithms  in  a  variety  of  environments 
where  precise  modeling  assumptions  are  violated.  Following  the  first  encouraging  results,  one  of 
our  objectives  is  to  develop  a  methodology  to  design  algorithms  which  are  capable  of  using  the 
available  a  priori  information  about  systems  to  be  identified,  estimated  and/or  controlled,  and  to 
avoid  limitations  imposed  by  preselected  parameterizations.  We  instead  plan  to  use  existing  physi¬ 
cal  parameters  and  thus  enhance  not  only  robustness,  but  also  serviceability  of  the  control  instru¬ 
mentation. 


SUMMARY  OF  RESEARCH: 


The  frequency-weighted  projective  controls  methodology  will  be  extended  to  include  distur¬ 
bance  rejection  as  well  robustness  to  unmodeled  dynamics.  This  will  involve  a  parameterization 
ol  classes  ol  dynamic  controllers  obtained  by  the  projective  controls  technique,  and  the  relating  ol 
this  parameterization  to  trequency  weightings  in  a  perlormance  index  The  parameterization 
should  summarize  the  design  freedom  available  at  a  given  stage  in  the  design,  and  provide  insight 
into  the  available  tradeoffs  of  perlormance.  robustness,  and  disturbance  rejection  Decentralized 
digital  control  will  continue  to  be  a  focus  in  the  problems  studied  Some  attention  will  be  given  to 
related  systems  modeling  issues,  such  as  choice  of  input-output  pairs  lor  decentralized  control  and 
the  obtaining  ol  corresponding  reduced-order  design  models  The  main  challenge  in  this  work  will 
be  to  formulate  in  a  usetul  way  the  controller  parameien/ations.  <:P.d  to  relate  these  cleverlv  to 
corresponding  frequency-dependent  weights  in  linear-quadratic  optimization  problems.  An  under¬ 
standing  of  the  fundamental  properties  of  such  optimization  problems  will  be  a  major  goal  of  this 
research.  This  should  result  in  a  noteworthy  increase  in  our  ability  to  deal  realistically  with  the 
control  of  high  order  lightly  damped  imperfectly  known  systems  using  decentralized  measurement 
and  computation,  and  in  the  presence  of  disturbance  inputs. 

The  development  of  new',  physically  motivated,  parameterizalions  will  be  a  major  research 
direction  in  our  adaptive  control  research.  The  goal  here  is  to  attack  the  design  of  reduced  complex¬ 
ity  controllers.  Robustness  and  convergence  properties  of  reduced-  complexity  adaptive  systems 
are  expected  to  be  superior  to  those  attainable  with  general-purpose  parameterizations.  such  as 
ARM A-parameterization.  Our  preliminary  results  summarized  in  the  proposed  research  section 
point  in  this  direction. 
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SCIENTIFIC  OBJECTIVE: 

Attempts  to  harness  several  newly  emerging  technologies,  which  lead  naturallv  to  large  scale 
computer  controlled  systems,  have  given  rise  to  a  significant  inieresi  in  t  he  underiving  field  ol  dis¬ 
tributed  and  decentralized  svstems.  Most  of  these  svstems  incorporate  task  decomposition,  mul¬ 
tilevel  coordination  and  control,  and  distributed  estimaiion  and  computation  Some  ol  them 
diverge  trom  classical  systems  in  that  the  basic  phenomena  governing  them  are  event-driven  and 
discrete  Hxamples  ol  such  event-driven  systems  are  many  and  include  automated  VLSI  produc¬ 
tion  systems  for  making  printed  circuit  boards  or  fabricating  walers  in  large  .illumes  in  flexible 
electronic  assembly  systems,  communication  networks  for  data  transmission  and  interconnected 
computer  systems.  They  require  descriptions  in  terms  ol  discrete  quantities  such  as  the  number  of 
PC  boards  to  be  retested,  or  other  dynamic,  logical  and  linguistic  .ariables  and  they  operate  as 
event-driven,  asynchronous,  largely  nondeterministic  processes,  and  hence  must  be  analyzed  as 
such  The  basic  goal  in  studying  such  large,  complex,  and.  typically,  distributed  and  decentralized 
systems  is.  firstly,  to  design  them  for  efficiency,  and  then  to  control  and  coordinate  them  to  attain 
this  efficiency.  Important  issues  that  arise  in  the  modeling,  control  and  coordination  of  such  svs¬ 
tems  are  information  flow,  learning,  aggregation,  time  scale  separation,  decentralized  estimation, 
and  distributed  computation.  The  prmctpal  objective  ol  this  research  unit  is  to  studv  such  issues  in 
order  to  enhance  our  understanding  of  the  behavior  and  control  of  large  scale  computer  controlled 
systems  in  a  distributed  network,  under  uncertainly  and  decentralization 


SUMMARY  OF  RESEARCH: 

During  the  first  phase  of  this  project,  one  of  our  locus  points  has  been  the  derivation  of 
optimum  control  laws  tor  distributed  systems  under  deceniraii/ed  none  lassie. il  inf  ormation  Such 
systems  and  the  ensuing  models  arise  in  practice  when  the  controller  memory  is  limited  requiring 
each  action  to  depend  on  the  most  recent  datum  Yet  another  scenario  would  be  the  situation  when 
the  physical  system  is  of  large  size  and  involves  various  subsystems  which  are  responsible  ndro 
dually  tor  specific  control  actions,  and  there  is  delay  in  the  transmission  ol  measurement  anil  .on 
trol  information  from  one  subsystem  to  another  In  this  context  we  have  shown  t ha l  a  sufx  iass  q 
such  problems  can  be  viewed  as  the  transmission  of  a  garbled  .ersion  ol  a  sequence  of  Gaussian 
random  variables  over  a  number  ol  noisy  channels  under  a  fidelity  criterion  furthermore  cr 
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have  shown  that  the  optimum  solution  (under  a  quadratic  loss  functional)  consists  of  linear! v 
transforming  the  garbled  message  to  a  certain  optimum  power  level  and  then  optimally  decoding  it 
by  using  a  linear  transformation  at  the  receiv  mg  end  We  have  also  extended  this  result  to  the  case 
when  the  communication  medium  has  unknown  statistical  description  and  the  optimum  policies  are 
derived  under  worst  case  performance  measures  involving  a  variety  of  fidelity  criteria. 

We  have  also  addressed  some  fundamental  questions  that  arise  in  the  control  of  discrete  event 
s>  stems.  One  particularly  important  issue  that  arises  in  modern  manufacturing  sv stems  is  the 
determination  oi  the  level  at  which  buffers  need  to  be  maintained  We  hav e  show  n  that  under  some 
circumstances,  a  run.  >-inveni<  <ry  policy  may  indeed  be  plimal.  even  in  the  lace  of  unreliabilu  \  Tins 
;s  a  surprising  result  since  buffers  are  usuallv  maintained  to  guard  against  unreliability  However 
there  is  ,i  deeper  phenomenon  at  work,  related  to  stead v  stale  inventor*  distnbuf ions  consisting  : 
bet:;  prohjbil;-;.  masses  as  well  as  absolute!'-,  continuous  components,  a  hich  maces  it  sometimes 
optimal  not  to  maintain  any  inventories. 
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SCIENTIFIC  OBJECTIVE: 

For  physical  devices  running  the  gamut  from  rotating  machinery  to  opto-electronic  laser  posi¬ 
tioning  systems,  nonlinearities  constitute  an  essential  part  of  any  reasonable  model.  Revolutionary 
developments  in  microelectronics,  coupled  with  recent  breakthroughs  in  the  understanding  of  the 
structural  properties  of  nonlinear  dynamic  systems,  open  up  new  possibilities  for  the  synthesis  and 
implementation  of  feedback  controllers  in  order  to  meet  the  steadily  rising  needs  for  high  perfor¬ 
mance  and  increased  efficiency.  However,  even  the  most  powerful  analytical  tools  that  are  available 
today,  such  as  exact  linearization  via  feedback,  nonlinear  input-output  characterizations,  and 
geometric  conditions  for  noninteracting  control  and  disturbance  rejection,  still  require  that  exact 
nonlinear  models  be  available.  This  is  in  direct  conflict  with  the  fact  that  the  mere  presence  of  non- 
linearities  makes  accurate  modeling  and  identification  difficult.  Consequently,  controllers  based 
upon  the  currently  available  theory  could  be  highly  unreliable. 

To  remedy  this  situation,  some  fundamentally  new  lines  of  research  are  needed  for  the  design 
of  realistic  robust  controllers  for  nonlinear  systems.  Driven  by  this  need,  our  research  program  has 
two  key  objectives.  The  first  is  to  develop  a  fundamentally  new  geometric-asymptotic  approach  for 
the  synthesis  of  robust  nonlinear  controllers;  this  approach  will  combine  the  advantages  of  exact 
geometric  notions  with  estimates  of  sensitivity  and  robustness  obtained  from  an  asymptotic 
analysis,  and  will  build  upon  synergisms  of  two  types  of  expertise  available  within  our  research 
'earn  The  second  objective  is  to  deal  with  circumstances  where  desired  levels  of  performance  can¬ 
not  be  achieved  by  feedback;  it  is  then  necessary  to  restructure  the  nonlinear  systems.  Our  objec¬ 
tive  here  is  to  develop  a  methodology  for  reconfiguring  nonlinear  systems  for  robustness  enhance¬ 
ment 


SUMMARY  OF  RESEARCH: 

In  the  past  year,  we  have  made  significant  inroads  in  the  study  of  robust  feedback  control  for 
nonlinear  systems.  We  have,  in  particular,  investigated  the  important  problems  of  the  robustness 
of  such  structural  properties  as  right-invertibility.  decouplabilitv  and  feedback  linearization  with 
respect  to  -.ampling  With  current  technology,  controllers  are  often  implemented  digitally  and  thus 
these  issues  are  among  the  most  basic  that  need  to  be  addressed.  In  this  connection  we  have  been 
able  to  show  that  right-invertibility  and  decouplabilitv  are  essentially  preserved  under  sampling 


but  that  feedback  linearizability  is  usually  destroyed  by  the  introduction  of  the  usual  sample  and 
hold  dev  ices  However,  by  passing  to  a  multirate  sampling  scheme,  feedback  linearizability  can  be 
recovered  Moreover  the  linearizing  feedback  can  be  computed  off-line  and  then  implemented  as  a 
table- look-up  scheme 

We  have  also  investigated  the  effects  of  possibly  nonlinear  modeling  uncertainty  on  a  variety 
of  control  system  synthesis  problems.  In  particular,  we  have  addressed  the  problem  of  robust  sta¬ 
bilization  ot  plant  models  exhibiting  structured  uncertainty  We  have  been  able  to  provide  a  new 
algorithm  called  "Residue  Iteration"  for  designing  optimal  controllers,  which  is  numerically  lar 
more  appealing  than  an.  existing  scheme  A  byproduct  of  this  research  is  a  new  algorithm  for 
ouier-t  unction  interpolation  that  v  leUts  much  lower  order  rational  mterpolants.  The  significance  of 
this  result  is  that  it  allows  the  sv  n thesis  ol  I ar  simpler  controllers  than  any  prev  ious  technique. 
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SCIENTIFIC  OBJECTIVES: 

Problems  involving  the  interaction  of  the  elements  of  a  multiple-user  communication  network 
are  among  the  most  important  and  most  challenging  problems  in  electronic  communications.  The 
performance  of  a  communication  network  depends  in  a  very  complex  way  on  the  routing  algo¬ 
rithm.  flow  control  mechanism,  acknowledgement  procedure,  channel  access  protocol,  error-control 
code,  signaling  scheme,  receiver  processing  method,  and  synchronization  technique  employed  in  the 
lower  three  layers  (network  layer,  data  link  layer,  and  physical  layer)  of  the  ISO  layered  model 
for  the  network.  The  objective  of  our  research  in  multiple-terminal  digital  communications  is  to 
gain  a  better  understanding  of  the  interplay  between  these  elements.  Our  research  will  focus  on 
the  issues  that  arise  in  mobile  radio  networks,  particularly  spread-spectrum  radio  networks,  which 
must  operate  in  hostile  environments  that  include  jamming  and  fading  This  imposes  additional 
requirements  on  the  network  in  terms  of  robustness  and  survivability  that  will  be  accounted  for  in 
our  research 

One  of  the  objectives  of  our  research  is  to  develop  new  signaling  methods  and  receiver  pro¬ 
cessing  techniques  that  will  exploit  the  implicit  and  explicit  redundance  that  is  present  in  the  sig¬ 
nals  and  messages.  Such  redundancy  exists  in  the  physical  layer  (diversity  transmission,  multipath 
signals,  and  modulation),  the  data  link  layer  (error-control  coding*,  and  the  network  lav er  (redun¬ 
dant  packets  and  messages).  We  are  particularly  interested  in  the  efficient  use  ot  implicit  and 
explicit  diversity  in  spread-spectrum  radio  transmissions  to  improve  «. ommumcation  pertormancr 
under  stressed  conditions  'e  g  jamming  or  heavy  network  traffic*  \nother  objective  is  to  develop 
network  protocols  that  are  compatible  with  and  take  advantage  ol  the  *  eat  ores  ol  spread  spectrum 
modulation.  Ot  particular  interest  are  algorithms  tor  distributed  scheduling  of  transmissions  A 
third  objective  is  to  examine  the  synchronization  problem  tor  spread-spectrum  radios  [  ffKient 
network  operation  requires  fast  acquisition  ol  the  spread  spectrum  signais  Moreo.er  the  acquisi 
lion  and  synchronization  svstems  must  be  able  to  operate  in  'he  presence  >1  multiple  transmissions 
jamming  and  lading  to  be  of  any  use  in  a  military  communications  net  work 
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SUMMARY  OF  RESEARCH: 

(  onsiderable  progress  was  made  on  the  problem  of  diversity  transmission  for  channels  with 
lading  and  burst  interference  Hurst  interference  can  be  either  pulsed  interference  or  partial-band 
interference  in  a  t requency-hop  (FH)  system  The  goal  of  our  research  on  this  topic  is  to  obtain 
diversity  combining  methods  that  perform  well  even  if  the  exact  nature  of  the  interference  and 
fading  are  not  precisely  known.  In  addition,  we  seek  combining  techniques  that  are  robust  with 
respect  to  the  received  signal  power,  the  interference  power,  and  the  thermal  noise  level 

In  [3)  we  have  demonstrated  that  the  ratio  statistic  provides  a  good  basis  for  diversity  com¬ 
bining  1  he  pert  or  ma  nee  of  a  rat  10- threshold  test  f  R  1  T  1  in  con  junction  with  both  sq  uare-!aw'  com¬ 
bining  and  majority-logic  combining  is  examined.  We  have  show  n  that  the  KTT  w  ith  majoritv 
logic  combining  gives  better  narrowband  interference  rejection  capahilin  than  the  RTT  with 
square-law  combining  for  diversity  of  order  three  or  greater  More  generally,  we  find  that  the  RTT 
is  an  effective  scheme  against  partial-band  interference  and  failing,  and.  when  used  in  conjuction 
w  ith  majority  logic  combining,  it  provides  good  performance  over  a  w  tde  range  of  channel  parame¬ 
ters 

In  other  research  on  HI  communication  in  the  presence  of  partial-band  interference  and  fad¬ 
ing.  we  have  begun  an  investigation  of  the  generation  and  use  of  side  information  with  Reed- 
Sclomon  and  errors- and-erasures  decoding.  Several  alternative  techniques  lor  obtaining  side  infor¬ 
mation  are  being  considered  and  research  on  the  reliability  of  candidate  techniques  is  in  progress. 
Generally  the  determination  of  the  reliability  of  the  side  information  and  the  effect  this  reliability 
has  on  system  performance  is  considerably  more  difficulty  for  channels  with  jamming  and  lading 
than  in  our  previous  work  on  side  information  for  channels  with  multiple-access  interference  onlv 
As  j  result  our  early  work  has  concentrated  on  narrowband  interference  rejection  where  some  of 
the  concepts  developed  tor  multiple-access  systems  apply  to  the  jamming  problem  as  well 

We  have  liegun  an  investigation  into  serial  techniques  for  the  coarse  acquisition  of  direct- 
sequence  spread-spectrum  signals  using  pseudo-noise  ( P\ )  sequences.  The  protocol  being  studied 
involves  correlating  the  received  signal  with  the  local  reference  signal  over  a  lime  period  (window) 
which  is  very  short  compared  to  the  period  of  the  PN  sequence.  A  three-way  decision  (Accept. 
Reject,  or  Expand  Window  Size)  is  then  made.  An  advantage  of  the  scheme  is  that  one  can  quite 
quicklv  decide  when  the  received  signal  and  local  reference  are  not  in  approximate  synchronism 
AV  e  have  used  a  new  bound  on  the  aperiodic  autocorrelation  of  PN  sequences  in  analyzing  the  per¬ 
formance  of  such  systems.  The  new  bound  is  much  lighter  than  previously  known  bounds  for 
sequences  id  long  period.  Our  results  indicate  that  the  average  time  to  achieve  coarse  synchroniza¬ 
tion  is  much  smaller  than  previously  determined.  However,  the  differences  are  most  apparent 
when  the  PN  sequences  have  very  long  periods. 

We  have  started  an  investigation  of  coding  for  channels  with  input  constraints,  such  as  dc- 
l  ree  coding  tor  optical  channels  and  magnetic  recording  subject  to  run-length  constraints.  We  have 
developed  an  encoding-decoding  algorithm  that  can  achieve  rates  up  to  channel  capacity  whether 
the  channel  is  noisv  or  not  In  the  noiseless  case  this  algorithm  is  easily  implementable.  requiring 
m.  low  precision  arithmetic  and  finite  memory  This  algorithm  can  also  be  used  for  dala- 
’ranslation  from  one  set  of  constraints  into  another  E  sing  this  algorithm,  we  have  been  able  to 
g  .e  ;n  -er'.arn  special  cases  a  constructive  proof  of  Orsiein  s  isomorphism  theorem  which  states 
•hai  anv  •  *0  Bernoulli  shifts  with  equal  entropy  are  in  a  certain  sense,  equivalent. 

AV  e  have  discovered  a  poiv  nomial-iime  algorithm  for  scheduling  transmissions  in  a  spread 
c  pec  ’rum  packet  radio  network  to  avoid  primary  collisions  of  packets  The  algorithm  produces 
ransmissu  n  schedules  which  maximize  the  rate  that  either  link  or  end-to-end  demand  is  satisfied 
n  '  he  oe'  work 
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SCIENTIFIC  OBJECTIVE: 

Statistical  signal  processing  functions  such  as  signal  detection,  estimation,  and  identification 
play  a  key  role  in  the  development  of  effective  communications,  radar  and  sonar  systems.  For 
example,  advanced  statistical  methods  are  emerging  as  being  particularly  important  in  digital  com¬ 
munications  systems  operating  in  channels  corrupted  by  interference  from  such  phenomena  as 
multiple-access  noise,  intentional  jamming,  and  impulsive  noise  sources.  Conventional  demodula¬ 
tion  methods,  such  as  coherent  matched  filtering,  often  suffer  serious  performance  degradation 
when  subjected  to  interference  of  these  types:  however,  this  degradation  can  frequently  be  elim¬ 
inated  through  the  use  of  more  sophisticated  signal  processing  techniques. 

A  central  issue  in  the  design  of  effective  signal  processing  procedures  for  systems  operating  in 
channels  such  as  those  noted  above  is  that  of  channel  identification.  Although  certain  aspects  of 
channel  identification  have  been  studied  extensively,  one  area  in  which  there  is  a  pressing  need  for 
further  research  is  that  of  identification  of  impulsive  channels.  Natural  impulsive  phenomena  are 
major  noise  sources  in  many  types  of  channels  including  ELF  electromagnetic  and  under-ice  acous¬ 
tic  channels.  Moreover,  man-made  impulsive  phenomena  are  a  principal  background  noise  source  in 
the  environments  in  which  military  radio  networks  must  operate.  Thus,  since  it  is  well-established 
that  impulsive  noise  can  be  extremely  detrimental  to  the  performance  of  communications,  radar, 
and  sonar  systems  if  not  properly  suppressed  and  since  impulsive  channels  often  exhibit  nonsta- 
lionarv  characteristics,  the  development  of  effective  techniques  for  identification  and  tracking  of 
the  characteristics  of  impulsive  noise  channels  is  an  important  problem  in  the  development  of  sys¬ 
tems  that  can  approach  the  performance  limits  set  by  such  channels. 

The  overall  objective  of  this  research  project  is  to  study  the  general  problems  of  identification 
and  tracking  of  impulsive  channels.  A  thorough  study  of  this  area  is  planned,  including  the 
development  of  suitable  channel  models,  the  derivation  and  analysis  of  optimum  batch  and  recur¬ 
sive  nonlinear  estimation  algorithms  for  identification/ tracking,  and  the  application  of  these  algo¬ 
rithms  to  develop  adaptive  techniques  for  the  reception  of  signals  passing  through  impulsive  chan¬ 
nels.  It  is  anticipated  that  the  results  of  this  study  will  find  application  in  a  broad  class  of  areas 
including  digital  communications,  sonar  and  radar. 
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SUMMARY  OF  RESEARCH: 

Our  research  efforts  during  this  reporting  period  have  focused  on  two  principal  aspects  of  the 
problem  ol  identifying  impulsive  channels.  These  are:  (1)  parameter  estimation  in  canonical 
statistical  physical  models  for  impulsive  channels;  and  (2)  probability  density  estimation  tor  gen¬ 
eral  parametric  channel  noise  models. 

Our  work  in  the  first  ot  these  areas  has  focused  on  the  Class  A  Middleton  model,  which  is  a 
physically  based  model  that  is  known  to  accuratelv  characterize  a  wide  variety  of  impulsive 
phenomena,  including  under-ice  acoustic  noise  and  man-made  electromagnetic  interl erence.  This 
model  has  two  basic  parameters.  A  and  T.  which  characterize  the  frequency  and  severitv  ol  the 
impulses  In  this  reporting  period,  we  have  explored  several  techniques  ter  estimating  these  param¬ 
eters  Irom  channel  measurements.  This  work,  w  hich  is  reported  in  [ !  ].  has  considered  batch  algo¬ 
rithms  only  and  thus  does  not  provide  methods  suitable  for  real-time  processing  or  for  tracking  ol 
time- varying  channels.  We  are  currently  exploring  recursive  versions  of  the  algorithms  studied  in 
[  1  ].  and  our  preliminary  simulation  analysis  has  shown  promising  results  for  these  new  algorithms. 

In  the  second  of  these  two  areas,  we  have  undertaken  a  theoretical  analysis  of  the  rates  of 
convergence  ol  probability  density  function  estimates  for  situations  in  which  the  density  function 
can  be  assumed  to  be  a  member  of  a  smooth  parametric  family.  These  convergence  rates  lead  to 
fieri ormance  approximations  that  are  applicable  to  a  wide  range  ol  stochastic  models  including,  in 
particular,  impulsive  noise  models.  Our  immediate  plans  in  this  area  include  the  application  of 
these  results  to  specific  parametric  models  for  impulsive  channels  (such  as  the  Middleton  Class  A). 
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SCIENTIFIC  OBJECTIVE: 

The  objective  of  this  research  is  to  tackle  fundamental  problems  of  electronic  materials,  dec 
ices  and  systems  in  a  timely  manner  and  to  provide  early  funding  on  the  start-up  of  projects  which 
present  immediate  opportunities  of  high  scientific  merit. 


SUMMARY  OF  RESEARCH: 

During  the  last  year  the  current  CS1  Director  ( W  k  Jenkins)  and  the  previous  ('SI  Director 
(T  V  1  rick,  presently  the  Head  of  Electrical  and  Computer  Engineering;  hav e  committed  a  total  of 
S250k  of  JSEP  discretionary  funds  toward  the  establishment  of  an  exciting  new  MBF  "Octopus- 
facility  that  will  be  housed  in  the  Coordinated  Science  l  aboratory  The  JSFP  commitment  is  in  the 
amount  of  Sl(X)k  for  the  first  year.  SI 00k  for  the  second,  and  S50k  tor  the  third  year  of  the 
current  JSFP  contract  A  team  of  researchers  from  the  University  of  Illinois  at  l  rbana-Champaign 
has  designed  a  unique  facility  in  which  state  of  the  art  surfaces,  interfaces,  and  multilavers  can  fie 
synthesized.  This  new  facility  will  lie  called  the  EpiCenter  an  abbreviation  for  the  L  ol  I  CSl 
Microelectronics  Joint  Center  for  Epitaxial  Growth  and  Surface  Science  The  S6M  initial  installa¬ 
tion  is  being  established  in  collaboration  with  the  Coordinated  Science  l  aboratory  and  the 
Microelectronics  Center  on  the  University  of  Illinois  at  l  rbana-Champaign  campus  As  detailed 
below,  funds  from  a  variety  of  sources  have  been  employed  to  purchase  the  equipment  now  on 
order. 

Seven  MBE  machines  in  the  EpiCenter  will  be  interconnected  by  uhv  t  5x10'“  torr )  stainless 
steel  transfer  lines.  From  the  same  vacuum  environment,  u  will  also  be  possible  to  access  advanced 
instrumentation  for  surface  modification  and  sample  characterization  The  crystal  growth  equip¬ 
ment  forms  a  cross  at  one  end  of  the  SO  x  40'  facility  Access  to  the  clean  environment  of  the 
facility  from  an  observation  area  is  attained  through  an  air  lock,  as  indicated  in  Figure  1  A  sam¬ 
ple  preparation  room  providing  class  lOO  clean  space  in  laminar  flow  hoods  is  accessible  onlv  from 
the  main  area.  A  lounge  is  provided  tor  planning  discussion  computation,  etc  during  facilitv- 
related  acliv  mes. 

Both  the  MBE  machines  and  the  uhv  transfer  lines  are  current  I  v  being  fabricated  bv  the  Per 
kin  FImer  Co  Figure  lb  shows  the  configuration  chosen  fur  the  v.  stem  Machines  1  and  2  m  Arm 
1  are  IU-\  compound  machines  with  2  adapted  to  gas-source  research  Morkot  Machines  1  and  4 

in  Arm  II  are  designed  for  xork  with  metals  1  3>  auh  both  eff  usion  and  e  beam  sources  and  1  <t 
polar  ceramic  materials  4i  Machine  5  is  contributed  bv  the  Lnr.ersit.  ol  lil.nuis  ai  l  roanu 
Champaign  Engineering  Research  Center  for  device  work  and  he  oher  machine  »>  n  Em  III  .> 
custom  designed  w  ith  e-beam  hearths  tor  research  mo  grow  th  d  thegruiplA  maier  als  Si  <  >e  ev 
Finally  ,  machine  *  on  Arm  !\  is  speciallv  adopted  to  Greene  s  research  using  on  beam  mixing 
select  particular  growth  paths  if  metasuble  structures  The  complex  hus  contains  a  x  -de  -arr  . 
it  crystal  growth  en  •  ironments  that  are  not  irdinaniv  .  >mpat  bte  x  if.  me  and  he:  t  -is  ill 
and  V  elements,  tor  example  are  electricallv  active  dopants  ,n  group  I'-  semiconduc'  rs  and  must 
therefore  be  emploced  with  selective  care  :n  group  I\  s  .  othesis 


U  special  feature  of  the  complex  is  that  the  seven  machines  are  irter.ernev.rei)  b-.  .it-.  Tjn.|v-r 
lines  held  at  a  pressure  of  5  x  10-11  torr.  The  purpose  is  to  permit  'rangier  t  a  gi.«*n  .uK-rjv 
trom  one  gross  th  chamber  to  a  second  without  surtace  contamination  1  hs  makes  u  r>.  sm(<  - 
gross  on  a  gi\ en  fresh  substrate  an  epitaxial  laser  sshich  is  m  t  v.'mpat Uve  *  :t rs  the  en .  lronme-o 
needed  to  'gross  the  substrate.  The  gross  th  ol  (iaAs  on  Si  is  a  particular  example  >:  an  impor.i"' 
ss  n thesis  problem  of  this  type  The  seven  machines  olfer  a  unique  tlesibilii  .  It  rehear. .m  •  ■!  rt;.- 
tspe  into  heteroepitaxs  In  particular,  the  Mill  and  sample  transier  eq  torment  t‘ir  agh« 
adapted  to  3"  ssalers  so  that  any  sample  mas  receise  success  is  e  pr.sexvinc  .n  .in  .  desired  >eu 
!  the  sesen  gross  th  chambers. 

3  s  anetv  ol  characterization  and  processing  instrument'  are  deigned  r.iO’  n  «»,.•  -rd  • 
t  tie  arms  These,  also,  are  lulls  adapted  tv'  the  3"  ssaler  size  »  ;  :  tie  x\  nttiesi'  .■  mpiex  1;;: 
mil'  in-sin .  modification  ol  samples  and  structure  charaderi/alun  xitmn  the  ^ 'motor.  .  a...  ..m 
ens  ironment  so  that  surface  contamination  is  avoided.  A  substrate  preparation  chamber  i'  .u'.whed 
to  the  end  of  Arm  1  so  that  substrates  can  be  o.xsdi/ed.  sputler-'d  annealed,  tested  eU  Opposite 
n  the  end  of  Arm  111.  a  rotating  anode  x-ray  source  and  ditf  ractometer  are  being  set  up  bs  /jhel 
tor  in  situ  studies  of  surfaces  and  multilayers  On  Arm  iV  a  I'M!  \PS  unit  is  being  installed  t • ' r 
surface  studies;  its  diagnostic  capabilities  are  enhanced  by  added  ion  scattering  spectroscope  ,  ultra 
Met  photoemission,  and  I.Ffil)  capabilities.  These  are  in  addition  to  the  RHI  I  I)  anv!  mass  --pec 
'r  cu'p'v  diagnostics  in  each  \1BT  machine.  Further.  MBH  machines  1  3  and  7  are  being  special1.  , 
-.-quipped  for  ellipsomelrv.  fluorescence,  and  other  optical  probes  ol  film  growth  I  flirts  are  nr* 
•vmg  made  to  obtain  Rutherford  backscatlering  equipment  for  interlace  and  surface  studies  n 
\rm  1\  of  the  s\  stem  A  planned  future  development  is  the  installation  of  focussed  ion  beam 
materials  modification  and  an  electron  microscope  on  Arm  1.  Plans  call  for  the  experimental  area  at 
•  r.e  ^r.vl  ■!  Arm  IV  to  become  a  surface  science  complex  to  which  I  reshiv  grow  n  surtaces  lrom  the 
.  mplex  .an  be  transferred  . 

\n  ut  Shddk  is  scheduled  from  University  of  Illinois  at  U  rbuna-('hampaign  DOF  Program 
.-vis  t  r  the  purchase  of  the  MBF  equipment  over  the  > ears  FY  S6  S7  and  SS.  The  remaining 
1  .-.vis  :  r  the  S6M  facility  have  been  contributed  by  \SF  Materials  Research  l  aboratory  t  S5»X»K  1 
\M  I  ngineering  Research  ('enter  (S200K),  JSFP  (S250K).  AFOSR  ( S6(K>K ).  and  the  U  niversity  of 
r  ‘is  j-  L  rbana-( 'hampaign.  together  with  a  gift  f  rom  the  Perkin  Timer  Corporation  of  uppru.xi- 
v  -  .  32  ' M 


w^lyjvlV^-lvlV. 


